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Abstract. We prove the existence of asymptotically flat solutions to the static vacuum Einstein 
equations on M — R 3 \ B with prescribed metric 7 and mean curvature H on dM ~ S 2 , provided 
H > and H has no critical points where the Gauss curvature K-, < 0. This gives a partial 
resolution of a conjecture of Bartnik on such static vacuum extensions. The existence and uniqueness 
of such extensions is closely related to Bartnik's definition of quasi-local mass. 



1. Introduction 



This paper is concerned with a conjecture of R. Bartnik [B3j . [B4] on the existence and uniqueness 
of static solutions to the vacuum Einstein equations with certain prescribed boundary data. To 
describe the conjecture more precisely, let M be a 3-manifold diffeomorphic to M 3 \ B where B is a 
3-ball, so that dM ~ S 2 . The static vacuum Einstein equations are the equations for a pair (g,u) 
consisting of a smooth Riemannian metric g on M and a positive potential function u : M — > 1R + 
given by 

(1.1) uRicg = D 2 u, Au = 0, 

where the Hessian D 2 and Laplacian A = trD 2 are taken with respect to g. The equations (jl.ip 
are equivalent to the statement that the 4-dimensional metric 

(1.2) g N = ±u 2 dt 2 + g, 
on the 4-manifold N = R x M is Ricci-flat, i.e. 

(1.3) Ric 9N = 0. 

This holds for either choice of sign in (jl.2p and since most of the analysis of the paper concerns 
the Riemannian data (g,u) in (jl.ip . we will assume g^ is Riemannian, and moreover identify t in 
(|1.2|) periodically, to obtain a metric on N = S 1 x M with t replaced by the angular variable 6. 

Given (M,g,u) as above, let 7 be the Riemannian metric induced on S 2 = dM and let H be 
the mean curvature of dM C (M, g) , (with respect to the inward unit normal into M) . Then the 
Bartnik conjecture [B4j states that, given an arbitrary such pair in C°°, 

(1.4) {~f,H) e Met°°(S 2 ) x C™(S 2 ), H > 0, 

there exists a unique asymptotically flat solution (g, u) to the static vacuum Einstein equations 
(jl.ip inducing the boundary data (7, H) on dM. 

This conjecture is a natural outgrowth of Bartnik's concept of quasi-local mass rngfO), [B2J, 
jB3j , defined as follows. Let (f2, 5) be a smooth compact 3-manifold with smooth boundary of non- 
negative scalar curvature, and define an admissible extension of (O, g) to be a complete, asymptot- 
ically flat 3-manifold (M,g) of non-negative scalar curvature in which (Q, g) embeds isometrically 
and is not enclosed by any compact minimal surfaces (horizons). Then 

(1.5) meffi) = inf{mADM(M) : (M,g) is an admissible extension of (Q,g)}, 
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where rriADM{M) is the ADM mass of (M,g), cf. [BlJ. Conjecturally, an extension (M,g) realizing 
the infimum in (|1.5p is a solution to the static vacuum Einstein equations (jl.ip on M = M \ £1 
which is Lipschitz, (but not smooth), across the junction d£l and for which the induced metric and 
mean curvature at the boundary of the interior and exterior regions agree: 

g\dM = g\dQ, Hqm = Hgti, 

leading to the boundary data (|1.4p . Observe that the boundary data (7,-ff) have the character of 
a mixed Dirichlet-Neumann type boundary value problem for the static equations (jl.ip . but the 
potential function u is absent from the boundary data. We point out that more standard Dirichlet 
or Neumann boundary data are not suitable for the (static) Einstein equations, cf. [An3j. 

In this paper, we prove the existence part of the Bartnik conjecture at least for a large class 
of boundary data (|1.4p . Thus, let Met m ' a (dM) be the space of C m ' a metrics on dM ~ S 2 and 
C r ^~ 1 ' a (dM) be the space of C m ~ l,a positive functions on dM. Let K 1 be the Gauss curvature of 
a metric 7 £ Met m > a (dM). 

Theorem 1.1. For any pair (j,H) G Met m,a (dM) x C™~ 1,a (dM), m > 3, satisfying the non- 
degeneracy condition 

(1.6) dH ^ on {K~ f < 0}, 

there exists an asymptotically flat solution (M, g) to the static vacuum Einstein equations (II. ip on 
M ~ M 3 \ B, realizing the data (7, H) at dM. 

The condition (jl.6p means that the mean curvature H has no critical points where < 0. In 
particular it holds for all metrics 7 of positive Gauss curvature on dM. 

To describe the approach taken towards the proof of Theorem 1.1, let £s = £™' a be the moduli 
space of AF static vacuum solutions (M,g,u) on a given 3-manifold M which are C m ' a up to dM, 
m > 3, and let 

(1.7) IL3 : £™ ' a -» Met m ' a (dM) x C^^idM), 

U B (g) = ( 7 ,H), 

be the Bartnik boundary map. We prove in Theorem 3.6 below that the space £s is a smooth 
(infinite dimensional) Banach manifold, and Hb is smooth and Fredholm, of Fredholm index 0. 
Let £ + be the open submanifold of £s for which the mean curvature H is positive, i.e. 

£+ = (Il B )- l {Met m > a {dM) x C+" 1,a (aM)). 

The Bartnik conjecture above may thus be rephrased to state that the smooth map LTs, restricted 
to £+, 

(1.8) II B : £ + -> Met m ' a (dM) x C^ a {dM), 

is surjective and injective, i.e. essentially a smooth diffeomorphism. 

A basic issue in the proof of the Bartnik conjecture is whether the map TIb on £ + is proper. As 
shown below, this is not true in general. However, let B nd C Met m,a (dM) x C T ^~ 1 ' a (dM) be the 
open subset on which the non-degeneracy condition (]1.6I) holds. Also let 

£ nd = (Tl B )~\B nd ). 

It is not known if the spaces £ + , £ nd or B nd are connected, so throughout the paper we assume that 
£+, £ nd and B nd are the connected components of these Banach manifolds containing the standard 
flat solution (M,gji a t, 1) where gji a t is the flat metric on R 3 \ -B 3 (l) with u = 1. The boundary 
data of this solution are thus (7+1, 2), where 7+1 is the round metric on the 2-sphere of radius 1. 



One of the main results of the paper, cf. Corollary 5.3, is then that the map TIb restricted to 

grid j g p r0 p er) [ e 

(1.9) U B : £ nd -> B nd , 

is a smooth proper Fredholm map. Such Fredholm maps have a mod 2 degree, degz 2 Ti-B, cf. [Sm|. 
Using the black hole uniqueness theorem [IJ , [RJ , [BMJ , we then prove in Theorem 6.1 that 

degi 2 TV B = 1, 

on £ nd . This proves Theorem 1.1, since a map of non-zero degree is necessarily surjective. 

We point out that the map TIb in (|1.8f) on the full space £ + is not proper. For example, for 
the flat solution gfi at with u = 1, there are boundaries given by embedded spheres (S 2 , 7^, ifj) 
with Hi uniformly positive, which converge smoothly in M 3 to a limit which is an immersed but 
not embedded sphere. Such a limit is then in the boundary d£ + , but the limit boundary data 
(7, if) £ Met m ' a (dM) x C™ _1,a ((9M); cf. Remark 4.3 for further discussion. This implies of course 
that the full boundary map n# in (jl.8p cannot be a diffeomorphism. 

In Theorem 5.1, the properness of the boundary map TIb is established for boundaries dM C 
(M, g) which are outer-minimizing, i.e. for which 

(1.10) area(S) > area{dM), 

for any surface S C M homologous to dM. (In fact, properness holds for solutions for which (jl.lOp 
holds in an arbitrarily small but fixed neighborhood of dM). However, (jl.lOp is not a property of 
the boundary data (7, if), since it involves the (global) structure of the solution (M,g,u). 

To tackle this problem, in Theorem 4.1 we prove that any static vacuum solution (M, g, u) € £ nd 
has a global foliation by surfaces St ~ S 2 with positive mean curvature and with So = dM. This 
is shown by proving that such solutions admit (a small perturbation of) a global inverse mean 
curvature flow (IMCF) starting at dM; this is where the non-degeneracy condition (|1.6j) is needed. 
Manifolds admitting such a positive mean curvature foliation automatically satisfy (jl.lOj) . 

Note that some condition, such as (II. 6j) . is needed to establish the outer-minimizing property. 
If (jl.lOp held for all solutions in £ + , then TIb in (jl.8p would be proper, which we know is not the 
case. 

Finally, it is worth mentioning that besides its relevance to general relativity, Theorem 1.1 
appears to be the first global (or semi-global) existence result for solutions of the Einstein equations 
on domains with a compact boundary component, in dimensions > 4. 

The contents of the paper are briefly as follows. In §2, we present background information on 
the structure of static vacuum solutions and choices of gauge. Section 3 discusses elliptic boundary 
value problems for the Einstein equations and proves the basic structure theorems on the moduli 
space £s and the boundary map II^. In §4, we prove that solutions in £ nd have a global positive 
mean curvature foliation. In §5, we then prove the requisite apriori estimates and establish the 
properness of TIb on £ nd . Finally, §6 contains the computation of the degree of lis and closes with 
several related remarks and open questions. 

2. Background Discussion 

Let M be a 3-manifold with compact boundary dM, and with a single open end E. (All of the 
results of this section and of §3 hold in all dimensions, but for simplicity, we restrict to dimension 
3). Apriori, dM need not be connected, although this will be assumed later on. As following 
(JOD-dOD, we let N = S 1 x M. 

Let Mets(N) = Met™' a {N) be the space of complete C m ' a static metrics on N, i.e. metrics of 
the form (|l.ip . m > 2. One has 

(2.1) Met™> a (N) ~ Met m ' a (M) x C™' a (M), 
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where C™' a {M) is the space of positive C m ' a functions on M. The space Eg = Eg(iV) of static Ein- 
stein (Ricci-flat) metrics on iV is equivalent to the space of pairs = (gM,u) £ Met(M) x C+(M) 
satisfying (jl.ip or (jl.3p ,(The smoothness indices will be occasionally dropped when unimportant). 

Recall that a complete metric g G Met m ' a (E) on an end E is asymptotically flat if E is diffeo- 
morphic to M 3 \ B, where B is a 3-ball, and there exists a diffeomorphism F : M 3 \ B — > E 1 such 
that, in the chart i 7 , 

(2.2) = 6ij + 0(r-\ d k9ij = 0(r~ 2 ), d k d e9ij = 0(r~ 3 ), 

in the standard Euclidean coordinates on IR 3 . The static vacuum equations (11.11) are invariant 
under multiplication of the potential u by constants. Throughout the paper, we assume that u is 
normalized so that u — > 1 at infinity, and that u is asymptotically constant in the sense that 

(2.3) u = l + 0(r- 1 ), d k u = 0{r~ 2 ), d k d e u = 0(r~ 3 ). 

Thus the 4-metric g^ is asymptotic to the product S 1 x M 3 at infinity. 

It is proved in [An2j that ends of static vacuum solutions (M, g, u) are either asymptotically flat 
or parabolic, where parabolic is understood in the sense of potential theory; equivalently, parabolic 
ends have small volume growth in that the area of geodesic spheres grows slower than r 1+e , for any 
e > 0. Moreover, asymptotically flat ends are strongly asymptotically flat in that the metric and 
potential have asymptotic expansions of the form 

, n .s 2m. . „ 2m 

(2.4) 9ij = (1 + — )<% + •••, u = 1 - — + ••• , 

where the mass m may apriori be any value m € M. These two behaviors, asymptotically flat and 
parabolic, are radically different and there is no curve of asymptotic structures for static vacuum 
solutions which join them. The finer behavior of asymptotically flat ends are described by the mass 
parameter m and higher multipole moments, cf. [BSJ. 

Let g be a fixed asymptotically flat (background) metric in Eg; henceforth Eg will denote the 
space of asymptotically flat static vacuum Einstein solutions. The static Einstein equations are not 
elliptic, due to their invariance under diffeomorphisms, and for several reasons one needs to choose 
an elliptic gauge. To begin, we consider the Bianchi gauge, and define 

(2.5) $3: Met™> a {N) -> S m ~ 2 ' a {N), 

$g(g) = Ric g + 5* g pg(g), 

where /% is the Bianchi operator, /3g(g) = 5g(g) + \dtr-g(g). Also, (5*X)(A, B) = \{{V A X,B) + 
(S7bX, A)) and 5X = —tr(5*X) is the divergence. The operator §g is a C°° smooth map into the 
space S m ~ 2 ' a (N) of static symmetric bilinear forms on N. 

Using standard formulas for the linearization of the Ricci and scalar curvatures, cf. [Be| for 
instance, the linearization of at g = g € Eg is given by 

(2.6) L(h) = 2(DQg)(h) = D*Dh - 2R(h). 

Clearly, L is elliptic and formally self-adjoint. In §3 we will discuss boundary value problems for 
and L. 

Next, the asymptotic behavior in the asymptotically flat end E requires the introduction of 
suitable weighted function spaces. We will use the standard weighted Holder spaces, although one 
could equally well use weighted Sobolev spaces. Thus, define Met™' a (N) C Met™' a (N) to be the 
subspace of metrics which decay to Euclidean data at a rate r~ s at infinity; more precisely, the 
component functions gtj and u of gjy should satisfy 

9ij - 5i d G C™> a (R 3 \B), u-16 C™> a (R 3 \ B). 
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Here C7 1 consists of functions v such that 



\v\\crn = ^supr fc+,5 |V fc ?;| < oo, 

k=0 



while C™' a consists of functions such that 



v U = v cp + sup mm r (x ,r (y ™+«+ d J LJ < oo, 

5 x,v \x-y\ a 

cf. [Blj . [LP| . Throughout the following, we assume the decay rate 5 is fixed, and chosen to satisfy 



(2.7) \ < 6 < 1. 

It is well-known, cf. [Blj . [LP] , that the Laplacian on functions, and Laplace-type operators on 
tensors, as in (j2.6fl . are Fredholm when acting on these weighted Holder spaces. 
The map <E> in (12. 5|) clearly induces a smooth map 



(2.8) $ : Met™> a (N) ,S™~ 2 ' Q (iV). 

Observe that <? is Einstein if $>g(g) = and flg(g) = 0, so that g is in Bianchi gauge with respect 
to g. (Although &g is defined for all g G Met™ ,a {N), we will only consider it acting on g near g). 

Given g G E s , let Me<™' a (iV) C Met™' a (N) be the space of C m ' a smooth static AF Riemannian 
metrics on N, satisfying the Bianchi gauge constraint 

(2.9) f3g(g) = at dN. 
As above, 

$ : M et™' a (N) -> ,S™~ 2 ' Q (A0. 
Similarly let be the space of metrics g G Met™' a (N) satisfying &g(g) = 0, and let 

(2.10) E c C Z c 

be the subset of static Einstein metrics g = gjy, Ric g = in Zc- The following result justifies the 
use of the operator $ to study Ec*. 

Proposition 2.1. Any static metric g = g^ G Zc sufficiently close to g is necessarily Einstein, 
g G Ec. Moreover, this also holds infinitesimally in the following sense. Let n be an infinitesimal 
deformation of g G Zc, i.e. k € KerDQ. If (3g(g) = 0, (for example g = g), then 

(2.H) Pg(lt) = 0, 

and n is an infinitesimal Einstein deformation, i.e. the variation of g in the direction n preserves 
(fO]) to 1 st order. 

Proof: Since g G Zc, one has &(g) = 0, i.e. 

Ricg + 5* g pg(g) = 0. 

Applying the Bianchi operator (3 g and using the Bianchi identity (3 g {Ric g ) = gives 

(2-12) &(<£(/%(</))) = 0. 

Setting V = /%(<?), a standard Weitzenbock formula shows that 2f3 g 5*(V) = D*DV — Ric(V). 

Also, since g,g G Met™' a {N), V G X^ s ,a {M), where X\+$ ,a {M) is the space of vector fields 

whose components are in C^ s 1 ' a (M). When acting on vector fields V with V = on dM, as in 

(|2.9p . the operator D*D is positive, with trivial kernel. Namely, if W G C™ +l ^ ,a is in the kernel of 
D*D, then integrating by parts gives 

/ \DW\ 2 + / (W,V N W) = 0, 

JB(r) JS(r) 
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where B{r) = {x G M : dist{x,dM) < r}. (Since W = on dM, there is no boundary term at 
dM). Letting r — » oo, the boundary integral tends to and so DW = 0, which in turn implies 
W = 0. 

Since D*D is self-adjoint and Fredholm, it has a smallest positive eigenvalue bounded away from 
0. For g sufficiently close to g, \Ric(V)\ = 0(e) pointwise and Ric(V) G C^\ s 3 ' a (M), so that we 
may assume that 2f3 g 5*(V) is a positive operator on V . Hence, again since V = on dM, the only 
solution of (I2.12p is V = 0, which implies g G Ec. 

To prove the second statement, let gt = g + tn. Applying the Bianchi operator j3 gt to $>(gt) gives 

(2-13) 9 M9t) = P 9t S* gt (Pg(9t)). 

Taking the derivative with respect to t at t = 0, one has (/? St 3> (<?<))' = + /?<&'. Both terms here 
vanish since g G and k is formally tangent to Zq- Hence the variation of the right hand side 
of (|2.13p vanishes. Since /%(#) = 0, this gives (3 g 5*(/3g(K)) = 0. The equation (12. lip then follows 
exactly as following (|2.12|) . with V = /%(k). 

■ 

Let denote the space of static diffeomorphisms of ./V which equal the identity on 

ON. Thus such diffeomorphisms decay to the identity at the rate r~ s and are independent of the t 
or ^-variable in (jl.2p . The group T>\ acts freely and continuously on Met(N) and E5 by pullback 
and one has the following local slice theorem for this action; we refer to [An3] for the proof. 

Lemma 2.2. Given anyg G E^' a and g G Met 1 J l ' a (N) nearg, there exists a unique diffeomorphism 
4> G close to the identity, such that 

(2.14) /%(</>*<?) = 0. 
In particular, (ffg G Met™' a (N). 

Lemma 2.2 implies that if g G E™' a is a static Einstein metric near g, then g is isometric, by a 
diffeomorphism in to an Einstein metric in E^' a , so that E^ ,Q! is a slice for K™' a under 

the action of T> r [ l+l,a . 

To prove that the moduli space £ is a smooth Banach manifold, (cf. Theorem 3.6), it is important 
to have a gauge with choice of boundary data in which the Einstein equations form a self-adjoint 
elliptic boundary value problem. This is not the case for the operator $ and we are not aware of 
geometrically natural self-adjoint boundary conditions for <£. For this reason, we will also consider 
another natural gauge, namely the divergence-free gauge. 

To do this, instead of consider 

(2.15) <!(<?) = $~ g (g) = Ric g - S -g + 5* g 5~ g g, 

where s is the scalar curvature of g = gjsf. The linearization of $ at g = g G Eg is given by 

(2.16) L{h) = 2(D$g) g (h) = D*Dh - 2R(h) - (D 2 trh + Sdh g) + Atrh g. 
In analogy to (|2.9p . define 

(2.17) Mei^ a {M) = {g G MetJ' a (N) : Sgg = at ON}. 

Similarly, let Z^' a = <t _1 (0) n Met 1 ^ ) ' a (M) and Ed C Zd be the space of Einstein metrics in 
divergence-free gauge with respect to g G E5. 

It is easy to see that Proposition 2.1 and Lemma 2.2 hold in this divergence-free gauge in place 
of the previous Bianchi gauge, with essentially the same proof. Thus = Zjj and for g G E^>, 



(2.18) 



5gg = on N. 
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Moreover, the diffeomorphism group T>\ transforms one gauge choice uniquely to the other. For 
instance, suppose (3g(g) = 0. Then we claim there is a unique <j) £ £q™ +1 > a suc h that 

(2.19) 5~ g {<P*g) = 0. 

At the linearized level, with g = g, this amounts to finding a vector field V with V = on dN such 
that if Pgh = then 5g(h + 5*V) = 0. This equation is equivalent to the equation 55* V = \dtrh, 
which is uniquely solvable for V with V = on dN. The local result in (|2.19p then follows from 
the inverse function theorem. 

3. The Moduli Space 

In this section, we study boundary value problems for the elliptic operators $ and <3?, and use 
this to prove that the moduli space 8s of static vacuum solutions is a smooth Banach manifold for 
which the boundary map Ub is Fredholm, of Fredholm index 0, cf. Theorem 3.6. 

We begin with the Bianchi-gauged Einstein operator <I> in (|2.5p . i.e. 

<5>g{g)=Ric g + 5* g (3g{g). 

Let A denote the 2 nd fundamental form of dM in M, A(X,Y) = {V X N,Y), where N is the unit 
inward normal into M, X,Y tangent to dM. Similarly, let Hm = trA denote the mean curvature 
of dM in M. Throughout the paper W T will denote the restriction or the orthogonal projection 
of a tensor W to T(dN) or T(dM). 

Proposition 3.1. Near any given background solution g E ,a , the operator <3? = ^ in (|2.5p with 
boundary conditions: 

(3.1) ftg(g) = 0, g\ T (8M) = 7M, H M = h at dN, 
is an elliptic boundary value problem of Fredholm index 0. 

Here the induced metric 7m is in Met m ' a (dM) while the mean curvature Hm of dM in (M, gj^) 
is in C m ~~ 1,a (dM). Note that the potential u does not enter this boundary data and so is formally 
undetermined at dM. Also the static property implies that (3g{g) vanishes in the vertical direction, 
/%(<?) W = o. 

Proof: It suffices to prove that the leading order part of the linearized operators forms an elliptic 
system. Recall from (|2,6p that the linearization of $ at g = g is given by 

L{h) = 2(D§g)(h) = D*Dh - 2R(h). 

The leading order symbol of L = 2D& at £' is 

(3.2) a(L) = -\a 2 I, 

where / is the Q x Q identity matrix, with Q = (n(n + l)/2) + 1; Q is the sum of the dimension 
of the space of symmetric bilinear forms on R n , together with the extra vertical S 1 direction. Here 
n = 3 but we give the proof for general dimensions. For static metrics, all components of the 
metric are locally functions on R n , and all derivatives in the vertical S 1 direction are trivial. In the 
following, the subscript represents the direction normal to dM in M, (or dN in N), subscript 
1 denotes the vertical direction, tangent to S , while indices 2 through n represent the directions 
tangent to dM. Note that one has h\ a = 0, for all The positive roots of (|3.2p are where 

£' = (f ,f), with multiplicity Q at £ G T*(R n ). 

Writing £' = (£,£«), i = 2, . . . ,n, (as above £i = 0), the symbols of the leading order terms in 
the boundary operators are given by: 

—2izh$k — 2i ijhjk + i£,ktrh = 0, k > 2, 

i>2 
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-2izhoo — 2i 6c^ofc + iztrh = 0, 



k>2 

h T = (7') T , and (H M )' h = u. 

This gives n+ n ( n ~ 1 ^ + \ = Q boundary equations, as required. Ellipticity requires that the operator 
defined by the boundary symbols above has trivial kernel when z is set to the root Carrying 
this out then gives the system 

(3.3) 2\Z\h ok -2iJ2Zjhjk + itktrh = 0, k>2, 

(3.4) 2\^\h 00 -2iJ2^khok-\^\trh = 0, 

k>2 

(3.5) h n = <t>, h T = 0, (H M )' h = 0, 

where cj) is an undetermined function. 

Multiplying (|3.3j) by i^k and summing gives, via (|3.5|) . 

k>2 

Substituting (|3.4p on the term on the left above then gives 

m 2 h o - m\ 2 trh = 0. 

Since trh = /loo + 4>, it follows that (ft = 0. 

Next, to compute H' M , we first observe that in general 

(3.6) 2A' h = V N h + 2Aoh- 25*(h(N) T ) - 5*{h 00 N). 

This follows by differentiating the defining formula 2 A = Cn9, and using the identities 2N' h = 
-2h(N) T - h 00 N, C N h = V N h + 2A o h. Since H M = tr M A, H' h = (tr M )' h A - tr M A o h and so 

(3.7) 2(H M )' h = tr M (V N h - 26*(h(N) T ) - 5*(h 00 N)). 

Hence the symbol of 2(HM)'f l is given by '^2k>2^' z ^ l kk — 2i£fc/iofc)- Setting this to at the root 
z = gives 

(3.8) YMhkk + 2i&fco*) = 0. 



k>2 



Via (13. 5p . this gives — 2iY^,k>2^khok = 0> and substituting this in (13. 4|) and using the fact that 
4> = gives 

2|£|/ioo - ICI^qo = 0, 

so that /ioo = 0. It follows from (j3.3l) that /lofc = and hence h = 0. This proves ellipticity 
of the boundary value problem (|3.ip and the Fredholm property follows from the fact that the 
Laplace- type operator L is Fredholm on Met™' a , cf. [LP]. 

Finally, it is straightforward to verify that the boundary data (13. ID may be continuously deformed 
through elliptic boundary data to elliptic boundary data for which L is self-adjoint and so of index 
0. This is proved in [An3] in a slightly different setting and the proof carries over here with only 
minor change, and so we refer to |An3] for further details. The homotopy invariance of the index 
then completes the proof. 



As noted in §2, we are not aware of a geometrically natural self-adjoint elliptic boundary value 
problem for In particular, the boundary conditions (13. ID are not self-adjoint. This property is 
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important for the proof of Theorem 3.6, and for this reason, we turn to the operator $ in (|2.15|) 
with linearization atg = g given by L in (|2.16p . 

Regarding boundary conditions for L, for h E S™~ ,a (N), let h T = h\dN an d l^ T ]o be the 
projection of h T onto the space of forms trace-free with respect to 7 = 7jv- Similarly, H' h denotes 
here the linearization of the mean curvature H = Hn of dN C N. 

We then have: 

Lemma 3.2. The operator L with boundary conditions 

(3.9) Sh = 0, [h T ] = 0, H' h = 0, 

is a self-adjoint elliptic operator. Moreover, under the first two conditions Sh = and [h T ]o = 0, 
the operator L is self- adjoint exactly for the boundary condition H' h = 0. 

Proof: It is a rather long (and uninteresting) calculation to prove that the operator L with 
boundary data f|3.9|) forms an elliptic system; this has been verified by computer computation 
using Maple. More conceptually, instead we will make use of Proposition 3.1 to simplify the proof. 
First, recall, (ADNj . |Tr| . that ellipticity of a boundary value problem is equivalent to the existence 
of a uniform estimate 

(3.10) \\h\\c<»,- < C(\\L(h)\\cm-2,a + WB^Wcm-i,* + \\h\\ C o) , 

where Bj is the part of the boundary operator of order j, together with such an estimate for 
the adjoint operator. As seen below, the boundary value problem is self-adjoint, so it suffices to 
establish (I3TTTH) . 

First, it is simple to prove (|3.10p for L in place of L via a slight modification of the proof of 
Proposition 3.1. Namely, for the boundary condition [h T ]o = 0, we have h T = <pj on dN, (in place 
of (|3.5p ). Note also that (|3.3p - (l3.4p hold, but without the trh terms. The analog of (|3.3[) then gives 

and hence, via the analog of (|3.4p . 

If | 2 />oo = -If IV, 

so that /too + = 0. Next, via the condition H' h = 0, the analog of (|3.8|) becomes 

5^(|£|fck* + 2i&/iofc) = 0, 

k>l 

which gives 

n|f|0 = -2iJ2Zkh 0k = 2|f|</>. 

Since n > 3, this implies = 0, and so /loo = 0, hence /iofe = 0. It follows that h = 0, which proves 
ellipticity of L with the boundary conditions f|3.9[) . Thus, (|3.10p holds with L in place of L. 
Next, one has 



(3.11) L = L — (D 2 trh — Atrh g) - 55h g. 
Thus to prove (|3.10p . it suffices to prove 

(3.12) \\5h\\ cm -i,* < C(l\L(h)\\ Cm -2, a + WBj^Wcm-^a + ||fc||oo), 

(3.13) \\D 2 trh\\ cm -2, a < C{\\L(h)\\ cm -2, a + \\Bj{h)\\cm-i, a + \\h\\ c o). 



From ([2TT5]) - (p36l) and the Bianchi identity, (as in (|2TT3|) ). one has 5L(h) = 255* (5h) and the 
operator 55* is elliptic with respect to Dirichlet boundary conditions. Since the boundary data 5h 
in (|3.9p is included in the boundary operators Bj, this proves (|3.12p . 
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Using this and taking the trace of (|3.1ip shows that 

\\D 2 trh\\ Cm -2, a < C{\\L{h)\\ Cm -2, a + ||.B,-(/t)||c«-J>« + \\N N (trh)\\ Cm -2,c, + \\h\\ c o), 

so that it suffices to prove that the boundary conditions B cover NN(trh). For this, a simple 
computation using (|3.7p . (cf. also (|3.19p below), gives 

(3.14) N(trh) = 2H' h - 5((h(N)) T ) - (5h)(N) + O(h), 

where 0(h) is of order in h. Using the standard interpolation | \h\ |cm-i,c« < e\ \h\ \c™,<* + e~ 1 ||/i|| c <o 
shows that it suffices here and below only to consider terms with the leading number of derivatives 
of h. 

Now the Gauss equations at dN are |j4| 2 — H 2 + s 7jv = s 9N — 2Ric(N, N) and hence, 
(\A\ 2 - H 2 + s w )' ft = -2L(h)(N, N) + 2S*6(h)(N,N) + 0(h). 

One has s'^ (h T ) = —Atrh T + 58(h T ) + 0(h T ) and A' h , H' h only involve first order derivatives in h. 

Writing then h T = B$(h) + ^trQ^h^N , it follows that trg^h at dN is controlled by L(h), Bj(h), 
in that 

\\trdNh\\c™^ < C(\\h\\ Cm -i,a, + ||L(/i)|| C m-2,a + \\Bj(h)\\ Cm -j, a ), 

and hence 

H^^HC" 1 ' — C(\\h\\(Jm.-l,a + | \L(h) I \(Jm.-2,a + | | Bj (h) \ \ Qm- j,a ) , 

i.e. h T is controlled at dN by L(h) and Bj(h). Next, at dN, one has -(5h)(T) = S7 N h(N,T) + 
V £t h(ei,T), which then gives control as above on (V jyh)(N,T), and so control on V N(h(N) T ). In 
turn, this gives then control on 5qn(\7 N(h(N) T )), which modulo lower order (curvature) terms, 
equals N(5(h(N) T )). The iV-derivative of (ETCH) also holds and shows that control of N(5(h(N) T )) 
implies control of NN(trh), so that f)3. 13|) holds, provided N(H' h ) is controlled. But the Riccati 
equation gives N(H) = — \A\ 2 — Ric(N, N); taking the linearization of this in the direction h shows 
that N(H' h ) is indeed controlled by L(h) and the boundary conditions Bj. This completes the proof 
of ellipticity. 

Next, we prove the operator L with boundary conditions (|3.9p is self-adjoint. To begin, integrat- 
ing the terms in the expression (|2.16p for L by parts over ./V gives 

/ (D*D(h),k)+ f (V N h,k)= f (D*D(k),h)+ f (V N k,h), 

JN JdN JN JdN 

[ 55htrk+ [ (5h)(N)trk = [ (h,D 2 (trk))- [ h(N,dtrk), 

J N JdN J N JdN 

and 

/ (Atrh)trk - / N(trh)trk = / (Atrk)trh - / N(trk)trh. 
Jn JdN Jn JdN 

Here the boundary terms on S(r) all tend to as r — > oo, since the components of h and k are in 
C™' a and 5 > \ . It follows that 

(3.15) / (L{h),k)+ [ (B(h),k)= [ (L(k),h)+ [ (B(k),h), 
Jn JdN Jn JdN 

where 

(3.16) {B(k), h) = (V N k, h) + h(N, dtrk) - (5k)(N)trh - trhN(trk). 
Setting Z(k, h) = (B(k), h) — (B(h), k), we thus need to show that 

(3.17) f Z(h,k)=0, 

JdN 

when h, k satisfy the boundary conditions (|3.9p , 
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Thus suppose h and k both satisfy f|3.9[) . A simple calculation shows that (5k) (T) = is 
equivalent to 

(3.18) (V N k)(N) T = 5 dN (k T ) - a(k(N)), 

where a(k(N)) = [A(k(N)) + Hk(N) T ], (all taken on dN), while (5k)(N) = is equivalent to 

(3.19) iV(fcoo) = 5 dN (k(N) T ) + (A, ft) - k 00 H. 
The same equations hold for h, and one also has 

(3.20) h T = (phi: an d k T = (f) k -y. 
We thus need to calculate 

B(k, h) = (V N k, h) + h(N, dtrk) - trhN(trk), 

and skew-symmetrize. To begin, write (Vjvfc,^) = ((V Nk)(N), h(N)) + {(V Nk)(ei), h(ei)), so that 
(V N k,h) = N(k 00 )h 00 + fa[N(trk) - N(k o)] + ^((^Nk)(N),h(N) T ), where we have used the 
relation tr^S/ Nk = tr^S/j^k — N(koo). Thus, B(k, h) equals 

(3.21) N(k 00 )h 00 +fa[N(trk)-N(k 00 )]+2((V N k)(N) 1 h(N) T )-N(trk)[^^ 
By (I3TT81) and (13301) . 

2((Vivfe)(iV), /i(iV) T ) = -2(# fc , ft(iV) T } - 2a(k, h) = -2fa5 dN (h(N) T ) - 2a(k, h) + u x 

where uji is a divergence term and a(k,h) = (a(k(N)),h(N) T ). Similarly, by $3l9\\ and <^M), 

N(k 00 ) = 5 dN (k(N) T ) + (fa - k 00 )H, 

where here and in the following 5 = 5qn. Note also that {h(N) T , dtrk) = trk5Q]y(h(N) T ) + UJ2, 
where UJ2 is another divergence term. Since (|3.17p involves integration over dN, in the following 
we ignore the divergence terms. Substituting these computations in (|3.2ip gives 

5(k(N) T )[h 00 - fa] + 5(h(N) T )[trk - 2fa] - (n - l)faN(trk) + H(fa - k 00 )(h 00 - fa) - 2a(h, k). 

When skew-symmetrizing, the last two terms H(<p k — fcooX^oo — (fth) — 2a(h,k) cancel, while the 
first three terms combine to give 

-(n - l)[fa5(k(N) T ) - fa8(h(N) T )] - (n - l)[N(trk)fa - N(trh)fa], 

or equivalently, (after dividing by n — 1), 

(3.22) - fa[N(trk) + <5(A;(iV) T )] + fa[N(trh) + 5(h(N) T )]. 
On the other hand, by (J3l6|) or (I3T71) . 

2(H') k = tr[V N k - 25*(k(Nf) - 5*(k 00 N)] 
= N(trk) + 25(k(N) T ) - k m H - N(k 00 )- 

Substituting (|3.19p gives 

2(H') k = N(trk) + 5(k(N) T ) - Hfa, 

so (|3.22p becomes 

-fa[2(H% + Hfa] + fa[2(H') h + Hfa] = -2fa(H') k + 2fa(H') h . 
This vanishes exactly when H' k and H' h vanish. This completes the proof. 

■ 

The main step in the proof of the manifold theorem, (Theorem 3.6), is the following result. 
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Theorem 3.3. Suppose m(M,dM) = and m > 3. Then at any g £ Kg, the map & is a 

submersion, i.e. the derivative 

(3.23) (D<% : T~ g Met^ a (N) - T^ (]f) S^ a {N) 

is surjective and its kernel splits in TgMet^ a (N) . 

Proof: The operator L = 2D&g is elliptic in the interior, and the boundary data in Lemma 3.2 
give a self-adjoint elliptic boundary value problem. Let Sg' a (N) be the space of C m,a symmetric 
bilinear forms on N satisfying the boundary condition B(h) = from Lemma 3.2, i.e. 

B{h) = {5h,[h T ] ,(H') h } = (0,0,0). 

Clearly, S™' a {N) C S™' a (N), where S™' a {N) = Tg(Met™' a (N)). Throughout the following, we 
set g = g. The operator L, mapping 

S™> a (N) -> S™- 2 ' a (iV), 

L{h) = f, B(h) = at 8N, 

is then Fredholm, of Fredholm index 0. On S^' a (N), the image Im(L) is a closed subspace of the 
range S rn ~ 2,a (N), of finite codimension, and with codimension equal to dimension of the kernel K. 

If K = 0, then L maps S^' a (N) onto ' a (N), which proves the result. Thus suppose 
Then as in (|3.15p . by the self-adjointness, one has for any h £ Sg' a (N) and k G K, 



(L(h),k) = / (h, L(k)) = 0, 

N JN 

since the boundary terms vanish and L(k) = 0. Thus Im(L\gm,a^) = K^~, (where K is taken 

with respect to the L 2 inner product). To prove surjectivity on S^ ,a (N), it thus suffices to prove 
that for any k G K, there exists h G Sp' a (N) such that 

(3.24) / (L(h),k)?0. 

JN 

Suppose then (|3.24|) does not hold, so that 

(3.25) / (L(h), k) = 0, 

JN 

for all h £ Sp' a (N), i.e. for which 5h = on ON. Integrating by parts, it follows that 



(3.26) / (h t L(k))+ / Z(h,k)=0, 

JN JdN 

for Z(h,k) as following (|3.16p . As before, the boundary terms at infinity vanish, since 5 > |- 
Choosing h S S^' a (N) arbitrary of compact support in N, it follows from (|3.26p that 

(3.27) L(k) = 0, 

i.e. k is formally tangent to Z = $ _1 (0). Of course this is already known, since k E K. Moreover, 
one also has 

(3.28) 5k = on N. 

To see this, let h = 5*V, with V any vector field vanishing on ON. Since g is Einstein and so 
(Ric - %g)' s * v = 0, it follows from {33]) and {MI) that L{h) = 5*Y, where Y = 255* V. As in 
Lemma 2.2, the operator 55* is surjective, (in fact an isomorphism), on vector fields vanishing at 
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dN, so that Y may be arbitrarily prescribed. Moreover, h € S^' a (N) if and only if Y = at dN. 
Then IET25]) gives 

0= / (L(<rn&>= / <^y,*:>= / (Y,Sk) + [ k(Y,N)= [ (Y,Sk), 
Jn Jn Jn JdN JN 

since Y = on dN. Here we have used again the fact that the boundary term at infinity vanishes, 
since \k\ = 0{r- 5 ) and \Y\ = 0{r- 1 - 5 ). Since Y is otherwise arbitrary, this gives (|3.28[) . 
Returning now to (|3.26p , (|3.27p gives 

(3.29) / Z(h,k) = 0, 

JdN 

for all h with 5h = on dN. Next, we choose certain test forms h € Sd(N) in (|3.29[) . First, choose h 
such that h = on dN. Then V^h is freely specifiable, subject to the divergence constraint 5h = 0; 
all computations here and below are at dN. Since h = 0, this constraint gives (V ]yh)(N) = 0, 
which is equivalent to the tangential and normal constraints: 

(3.30) (V N h)(N,T)=0, 

(3.31) N(h 00 ) = 0, 

for any T tangent to dN. Choosing h and Vtv^ satisfying h = and (|3.30j) - (|3.31j) at dN, the 
terms (Vat/i)(Ti, T2) are freely specifiable on dN, where T\,T<i are any vectors tangent to dN. 
Substituting such h in (|3.29p and using (|3,28p . it follows that 

(3.32) / (Vjv/i, k) + (k 00 - trk)N{trh) = 0. 

JdN 

Now choose Vjv/i = fg T , where g T = g\x{dN)- This choice satisfies the constraints (j3.3U|) - (|3.31|) . 
The integrand in (|3.32p then becomes ftr T k — N(trh)tr T k. Since N(trh) = (Vat/i,<7) = nf, and 
since / is arbitrary, it follows that tr T k = 0. In turn, since the tangential part of Vjyh is arbitrary, 
(HEMP implies 

(3.33) k T = 0, on dN. 
Lemma 3.4. At dN, one has 

(3.34) (A' k f = 0, 

i.e. (V N k) T = 2[5*(k(N) T )] T + k 0Q A, since k T = 0, cf. (IBTBj) . 

Proof: The proof is a straightforward, but rather long computation. To begin, as preced- 
ing (|3.2ip and using QESD , one has (V N h,k) = 2({V N h)(N), k(N) T ) + N{h 00 )k 00 . By (|348l) . 
(X7 N h)(N) T = 5 dN (h T ) - a(h(N)), so that 

(3.35) I (V N h,k}= [ 2(5 dN (h T ),k(N) T ) + N(h 00 )k 00 -2a(h,k) 
JdN JdN 

= [ 2(h T ,(5 dN y(k(N) T )) + N(h 00 )k 00 -2a(h,k). 
JdN 

Further, for Z tangent to dN, one has (5 9N )*{k(N) T ){Z, Z) = {V T z k{N) T ,Z) = (V z k(N) T ,Z) = 
5*(k(N) T )(Z, Z), where now 5* is taken with respect to the ambient metric g^, (not the boundary 
metric 7^). So this gives 

(3.36) f (V N h,k}= [ (h T ,25*(k(N) T )) + N(h 00 )k 00 -2a(h,k). 
JdN JdN 
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On the other hand, one computes (V N k,h) = ((V N k) T ,h T ) + (V N k(N),h(N)) = ((V N k) T ,h T )- 
(a(k(N)), h(N) T ) + N(k 00 )h 00 , again by (Engl) and (ET55J) . Taking the difference of this with (ET351) 
and noting that a is symmetric, gives 

(3.37) f (h T ,(V N k) T -25*(k(N) T )) + N(koo)h 00 -N(h 00 )k 00 = E, 

JdN 

where via (j3.16p - (j3.17p . E is given by 

E = / [fe(iV, dtr/i) - h(N, dtrk)] - [N(trh)trk - trhN(trk)]. 
JdN 

Computing this term-by-term gives: kooN(trh) + (k(N) T \dFtrh) — hooN(trk) — (h(N) T \dFtrk) — 
N (trh)trk + trhN (trk) . Since trk = koo, the first and second-to-last terms cancel. Integrating over 
dN and using the divergence theorem shows that 

(3.38) E = I trh5 9N (k(N) T ) - k 00 5 T (h(N) T ) - h 00 N(trk) + trhN(trk). 

JdN 

Next we claim that 

(3.39) 5 dN (h(N) T ) = N(h 00 ) + Hh 00 - (A, h), 

and similarly for k. This follows from the following computation: 5QN(h(N) T ) = 5g]y(h(N)) — 
5 T (h 00 N) =6 dN (h(N)) + Hh 00 , while S dN (h(N)) = 5(h(N)) + N(h 00 ). Since 5(h(N)) = (5h)(N)- 
(A, h), this gives the claim. Substituting (|3.39p into (|3.38p . and using {A, k) = implies that 

E= trh{N(k 00 ) + Hkoo) ~ k o(N(h 00 ) + Hh 00 - {A, h)) - h 00 N{trk) + trhN{trk), 
JdN 

and rearranging terms gives 

(3.40) E= (A,h}k 00 + N(trk)[trh-ho }+trhN(k 00 )-kooN(h 00 ) + H(trhk 00 -trkhoo). 

JdN 

Now substitute (|3.4U|) into (|3.37|) : the kooN(hoo) term cancels to give 

(3.41) f {h T ,(V N kf-25*(k(N)f)-(A,h)k 00 = 

JdN 

- / N(k o)h 00 - N(trk)[trh - h 00 ] - trhN(k 00 ) - H(trhk 00 - trkh 00 ). 
JdN 

The integrand on the right combines to: —N(koo)(hoQ — trh) — N(trk)[hoo — trh] — Htrk(hoo — trh) = 
— [N(koo) + N(trk) + Htrk](hoo — trh). Since hoo — trh = —tr T h = —(h T ,g T ), and since h T may 
be chosen arbitrarily, (the constraint 5h = imposes no constraint on h T ), it follows that 

(3.42) (VivA:) T = 2[5* (k(N) T )f + k 00 A + [N(k 00 ) + N(trk) + Htrk]g T . 
To complete the proof of (|3.34p . we thus need to show that 

(3.43) N(k 00 ) + N(trk) + Htrk = 0. 

To obtain (EjjSD , take the 5 T -trace of (ET4"2l) . One has (V N k,g T ) = N(trk) - N(k 00 ), while 
(5*(k(N) T ),g T ) = (Ve^AOV*) = (V e MN),ei)-k 00 H = ((V ei k)(N),e l )-k(V ei N,e i )-k 00 H = 
—N(koo) — kooH, the last equality using (|3.33|) and (|3.28|) . This gives 

N(trk) - N{k 00 ) = -2N(k 00 ) - 2k 00 H + k 00 H - n[N(k 00 ) + N(trk) + Htrk], 

which implies (|3.43|) . This completes the proof of the Lemma. 

■ 
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To complete the proof of Theorem 3.3, (|3.33|) and (|3.34[) show that 

k T = (A' k ) T = 0, 

at dN. One also has L(k) = 5k = on N, so that k is an infinitesimal Einstein deformation on N. 
By the local unique continuation result of AHJ, together with the global hypothesis vri(M, dM) = 0, 
it follows that k = 0. This shows that L is surjective. The fact that its kernel splits is standard, 
cf. [An3] . This completes the proof. 

■ 

Via the implicit function theorem, one obtains: 

Corollary 3.5. Suppose vri(M, dM) = and m > 3. Then the local spaces E^' a are infinite 
dimensional C°° Banach manifolds, with 

(3.44) TgE D = Ker(D$~). 

Proof: This is an immediate consequence of Theorem 3.3, the fact from Proposition 2.1 that 
Ex> = Zd, (cf. (|2.18|) ). and the implicit function theorem, (or regular value theorem), in Banach 
spaces. ■ 

This leads to the main result of this section. 

Theorem 3.6. Suppose 7Ti(M,dM) = and m > 3. Then the moduli space £s 
smooth infinite dimensional Banach manifold for which the boundary map 

(3.45) n B : £ S -> Met m ' a (dM) x C m ~ l ' a (dM), 
is a C°° smooth Fredholm map, of Fredholm index 0. 

Proof: Recall from §1 that the moduli space £g of static vacuum Einstein metrics is defined to be 
the quotient E™' a /T>™ +1,a . The local spaces E_d are smooth Banach manifolds and depend smoothly 
on the background metric g, since the divergence-free gauge condition (12.180 varies smoothly with 
g. As noted preceding Lemma 2.2, the action of T>\ on E is free and the local spaces Ex> are smooth 
local slices for the action of T>\ on Eg. Hence the global space Eg is a smooth Banach manifold, as 
is the quotient The local slices E_d represent local coordinate patches for £"5. 

Proposition 3.1 implies that the boundary map U B : E™'° -» Met m ' a (dM) x C m - l > a (dM) 
is smooth and Fredholm, of Fredholm index 0. Moreover, lie is invariant under the action of 
P™ +1 ' Q (M) on £™ ,a and so it descends to a smooth Fredholm map as in (|3.45p . still of index 0. 



£ 
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4. Positive mean curvature foliations. 

In contrast to the discussion in the previous section which held in all dimensions and for al- 
most arbitrary topologies, we assume here that dimM = 3 and further that dM ~ S 2 with M 
diffeomorphic to R 3 \ B, where B is a 3-ball. 

Recall from §1 that £ nd C £™ ' a is the space of static vacuum solutions (M, g, u) for which the 
boundary dM has positive mean curvature and satisfies the nondegeneracy condition (jl.6p . 

\VH\ ^ on {K < 0}. 

Throughout the rest of the paper, the connection on M will be denoted by V and the induced 
connection on a surface ScM will be denoted by V. 

The purpose of this section is to prove that solutions (M, g, u) G £ nd have outer-minimizing 
boundary (jl.lOp . i.e. 

(4.1) area(S) > area(dM), 
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for any surface £ enclosing dM. This is established by proving that all solutions in £ admit a 
foliation £j by surfaces of positive mean curvature with So = dM. It is well-known, (via a standard 
maximum principle argument), that the existence of such a foliation implies (|4.ip . 

To begin, for technical reasons it will be convenient to have a more precise lower bound for 
the Gauss curvature at critical points of the mean curvature. Let e : £ nd — > M + be a small error 
function, tending to sufficiently fast at the boundary of £ nd ; thus e < eq, where Eq is a fixed 
small constant, and for any £o > s > 0, the set of (M,g,u) for which e(M,g,u) > s is a 'weakly' 
compact set in £ nd , i.e. compact in the C m,a norm, a' < a. Now given a positive constant wq < 1, 
define S™* C £ nd to be the space of static vacuum solutions with boundary having positive mean 
curvature and the following additional properties at dM: 

(4.2) |Vlog#f + eK > Wq 1 -1 + Oe, and 

1 < w Q 1 < 1 + e, 

where 9 = ^min||y^| =0 j \K\ > 0. Observe that given any solution (M,g,u) G £ nd , there exists wo 
with Wq 1 — 1 sufficiently small, such that (M, u) G f . Thus £ nd = U Wo £™ d . It is clear that the 
spaces £™o are open submanifolds of £"5. 

Next we consider solutions (M, u) which admit foliations with leaves of positive mean curvature 
which are small perturbations of the inverse mean curvature flow (IMCF) starting at dM. Such 
foliations have leaves Ef ~ S 2 given by level sets of a smooth function t : M — > R, (it 7^ 0, which in 
the 0-shift gauge satisfy the evolution equation 

(4.3) d t = (wH^N. 

Define then £p d Wo C £™ d to be the subset of solutions which possess a smooth generalized inverse 
mean curvature flow (|4.3|) with w G C m (M), which at dM satisfy w\qm = wo, ^ 1 ^w\qm = for 
2 < I < m, V nw^^m = —eH~ l so that Au> -1 |,9m = — e, and with the following additional 
properties: on all of M 

(4.4) 1< iv' 1 < 1 + ee~\ \V l w~ x \ < ojee - *, Aw" 1 < -e 2 e _2i 

where {a^}^]^ is an appropriately chosen increasing sequence, and on each leaf £$ the nondegeneracy 
condition holds 

(4.5) \Vlog{wH)\ 2 + eK > uT 1 - 1 + ee ~*. 

It is clear that fp^Q is an open submanifold of £% d . As in §1, we assume here and throughout the 
following that we are taking the connected component of the standard solution (R 3 \i? 3 (l), gjiat, 1) 
in each of the relevant spaces above. The next result proves the statement (|4.ip . 

Theorem 4.1. The space £p d w is closed in £^, d , and hence 

(4-6) £ n F d WQ = £ n w d . 

Consequently, every solution (M, g, u) G £ nd admits a foliation E t of positive mean curvature, with 
E = dM. 

Proof: The strategy may be described as follows. Suppose that (M, gi, Uj) is a sequence of static 
solutions in £ F wo which converges to a limit (M, g, u) in £™ d in the C™' a topology. By assumption 
each member of the sequence has a foliation Ti of positive mean curvature, and a function Wi 
used to construct the foliation. According to the estimates (|4.4|) . the sequence u>i converges (in a 
subsequence) to a function w. Although the top order derivatives of w are not necessarily 

continuous, they are appropriately bounded, which is enough for all the arguments to follow; at 
the end, to regain the full regularity one may use a C m approximation. Moreover the structure 
of the foliations naturally leads to a positive lower and upper bound for the mean curvature in 
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terms of the ambient geometry. With this estimate and the bounds for Wi, the full geometry of the 
foliation may be controlled. Thus Ti will converge (in a subsequence) to a C m ~ 1,a limit foliation 
J- of (M,g,u). However the limit foliation may have points at which equality holds in (|4.4|) and 
(|4.5p . so the last step entails making a small perturbation to regain the strict inequalities. 

We begin by obtaining a uniform lower bound for the mean curvature using an energy method 
associated to the flow (|4.3p . The same method gives also a uniform upper bound for H. Let F(u) be 
an arbitrary smooth positive function of the potential, let f(t) > and tp(x) > also be arbitrary 
but smooth, (all three will be chosen below). We set Rnn = Ric(N,N). Then by a computation 
given in the Appendix, 



(4.7) 



f- x d t ( iP(x){f{t)wHF{u)) x = I XwiPF\wH) x - 2 Q, 

J-Et JT, t 



where 

(4.8) Q = -(\-l)\Vlog(wH) + (logF)'uVlogu\ 2 -((logF)" + (logF)' 2 )u 2 \Vlogu\ 2 
-\A\ 2 - (1 - (log F)'u)R NN + X^H 2 + HN(logw) + 2H (log F)'uN (log u) 
+w(log f)'H 2 - (Vlog(u>V), (Vlog(wH) + (logF)'iiVlogu)}. 

We will let A — > — oo, (and later A — » +oo), and will choose the unspecified functions appropriately 
so that the right side of (|4.7p is negative. The choice of the unspecified functions will vary depending 
on whether the region in question contains critical points of log(wH). In order to restrict the first 
variation formula (|4.7p to these different types of domains, we will employ ip as a cut-off function. 

To begin consider a region Q x (s,t) on which \V log(wH)\ ^ 0, and let ijj £ C^°(Q). Then by 
choosing F(u) = 1, all terms in Q will be dominated by —(A — 1)[V log(wH)\ 2 as A grows large, 
except perhaps the last term. However the last term may be estimated by 



(4.9) | (Vlog(io^),(V logoff) + (logF)'nVlogu))| 

< HAl^lVlog^VOI 2 + ||A||Vlog(^) + (logF)' 

Therefore, using the term w(log f)'H 2 in (|4.8p and by choosing 

wip(wH) x ^ 2 \Vlog(wip)\ 2 



(4.10) (log/)' = |A|- 



the first variation will be negative. 

Now consider a portion of the set {|Vlog(u;iT)| = 0}. The nondegeneracy condition (|4.5p implies 
that we may choose a domain x (s, t) containing this portion and on which K > ^6e~ t . As above 
we again choose ip £ C£°(£l) to restrict attention to this domain. Here we can no longer rely on 
\X7log(wH)\ to dominate terms in the first variation formula. Thus we must examine the remaining 
terms more closely. To this end, the Gauss equation on gives \A\ 2 — H 2 + 2K = —2Rnn, and 
hence one has 

-\A\ 2 -(1- (log F)'u)R NN = -1(1 + (log F)' 'u)\A\ 2 + (I- (log F)'u)K 
(4.11) -Ul- (log F)'u)H 



2 



This along with the coefficient of | V lo gn| 2 in gSD sug gests that one should choose F(u) to satisfy 
(4.12) (log F)" + (log F)' 2 < 0, 1 + (log F)'u < 0. 

Both of these inequalities will hold if (logF)' = (2u - c)~ l (so that (logF)" = -2(logF)' 2 ) with 
the constant c satisfying 2u < c < 3u. Notice that it is not necessarily possible to find a constant 
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having this property globally; however by restricting the size of O x (s, t) appropriately such a c 
will exist for the domain in question. The remaining terms in Q form a polynomial in H: 

(4.13) [wQog /)' - |(1 - (log F)'u) + A" 1 ] H 2 

+ [N(\og w) + 2(log F)'uN{\og u)\H 

(log F)'u)K. 

In light of the lower bound for the Gauss curvature, we may choose 

(4.14) (lo g/ y = C+ |A|-^ IJJ p)x 

for some appropriately large constant C > (depending only on w, u, and the Bartnik data), to 
guarantee that (|4.13|) is positive. Moreover the last term in (|4.14p is used as before to dominate 
(|4.9p . It follows that the first variation formula is again negative. 

We have shown how to construct the unspecified functions in domains $7 x (s,t), if f2 and t — s 
are sufficiently small. Therefore 

/ (f(t)wHF) x < f i,(x)(f(t)wHF) x 



and it follows that 



L 



< / 4>(x)(f(s)wHF) x 

< / (f(s)wHF) x , 



(f(»)<nIIF) x ) ( [ (f(t)wHF) x 



's s / \JT, t n{4>=i} 

Note that the cut-off function ip is chosen to take the value 1 on a compact set of f2, and in such 
a way that (|4.10p and the last term in (|4.14p converge to zero as A — > — oo. Thus upon taking the 
limit we obtain 

(4.15) mm f(s)wHF(u) < mm f(t)wHF(u), 

s a E t n{'0=i} 

where f(t) is either 1 or exp(Ct). 

The desired lower bound for the mean curvature may now be obtained as follows. Cover dM 
by appropriate domains O and construct the functions F, f, and ip as above according to whether 
or not Q contains critical points of log(wH). By applying (|4,15p we may then estimate mins t H 
in terms of mmg^ H for small t. Next repeat this procedure with dM replaced by S^. Hence by 
working outward we eventually find that for all t 

(4.16) mini? > C(t) > 0, 

St 

where the function C(t) depends only on the ambient geometry (M,g,u) and the function w. By 
letting A — * oo, the same arguments give a uniform upper bound maxs t H < K{t) < oo. For this, 
one changes the sign in (|4.12p . so that (logF)' = — (2u — c) -1 , with 2u < c < 3u, as well as the 
signs for (log/)' in (|4.10[) and (|4,14j) . 

Hence the lapse (wH)^ 1 of a foliation T as above is uniformly controlled by the ambient geometry 
(M,g,u) and the given bounds on w. This applies in particular to the sequence of foliations T% 
preceding (|4.7p . 

Moreover, with these bounds for the mean curvature in hand, it is straight forward to obtain a 
bound for the second fundamental form by following Smoczyk's arguments [SmoJ, (see also [HI]). 
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More precisely there exists a constant C depending only on the ambient geometry, the bounds for 
the mean curvature, minu>, |w;|c2, and maK.QM(wH) 2 (l + \A\ 2 ) such that (wH) 2 \A\ 2 < C(l + t 2 ). 
Furthermore by applying the estimates of Krylov [K] in the standard way, we may use the C° bound 
for the second fundamental to obtain control over all higher order derivatives: there exist constants 
C(l,f3,t) depending only on the ambient geometry, minu;, and such that l-Al^j-i.p < 

C(l,/3,t), I < m. Hence the foliations J-{ converge appropriately (in a subsequence) to a limit 
foliation T . 

As described earlier, the limit foliation T and function w may admit equalities in (|4.4p and (|4.5p . 
However here we will show that it is possible to perturb the foliation to regain strict inequalities. 
Each leaf E$ will be deformed in the normal direction by an amount /izt, where \i is a small 
parameter. For convenience we will drop the subscript t when there is no danger of confusion. The 
first task is to compute the new lapse function. Let v = (wH) -1 be the old lapse and let tilde 
signify a quantity after perturbation. Then up to first order the new lapse is given by 

v + tidtz = (uT 1 + iiHdtz)R- x 

= (uT 1 + ixHd t z){H~ x - iiDH-\z\) + 0(/i 2 ) 

= (w' 1 + ixHd t z - ixw^HDH- 1 ^})^- 1 + 0(/i 2 ). 

It follows from the formula for the variation of H, (2 nd variation of area), that 

(4.17) fiT 1 = w' 1 +i2{Hd t z- W - 1 HDH- 1 [z])+0{n 2 ) 

= w- 1 + fi[Hd t z - [wH)~ x (Kz + {\A\ 2 + R NN )z)] + 0{fi 2 ). 

In order to ensure that (|4.4p and (|4.5p regain strict inequality after perturbation, we choose z in 
the following way. Let x £ C°°(M) satisfy x = 0{e 2 e~ 2t ) (this function will be specified in (|4.2ip 
below), and let (j> G C°°(M) be a cut-off function agreeing with w~ l — x to all orders at dM and 
vanishing for t > to where to is small and depends on the ambient geometry. Then define z to be 
the solution of the Cauchy problem 

(4.18) Hd t z-(wH)- 1 (Az + (\A\ 2 + R NN )z) = <i> + x-w- 1 on M, z\ dM = 0. 

The choice of <p implies that z vanishes identically at the boundary, and hence u> _1 agrees there 
with w~ 1 to all orders. 

Lemma 4.2. Equation (I4.18f) admits a unique solution. Furthermore there exists a constant C 
depending only on the ambient geometry such that 

(4.19) sup(|*| + \Vz\ + |V 2 z| + \d t z\ + \Vd t z\) < Ce t/2 . 

s t 

Proof: When written in a more conventional form, equation (14.150 becomes 

d t z - w- x H-*Kz - w^H^dAf + R NN )z = H' 1 ^ + x - w' 1 ). 

This is a linear parabolic equation and thus the Cauchy problem admits a unique solution for all 
time. The bounds (14.191) will follow from the maximum principle and Schauder estimates. First 
observe that by setting z = (1 + tYe 1 ^ 2 ^ we obtain 

(4.20) d t z - w^H^Az - w^H' 2 (\A\ 2 + R NN - (± + <r(l + t)' 1 ) wH 2 ) z 

= (l + tJ-'e-'^ff-^ + x-uT 1 )- 

Under IMCF (or the current small perturbation of IMCF) the mean curvature and second funda- 
mental form fall-off like e _i//2 , and furthermore the leaves approximate round spheres at infinity. 
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Therefore if a > then the zeroth order coefficient satisfies 

\A\ 2 + R NN - (i + a(l + t)- 1 ) wH 2 
= -a(l + t)~W 2 + (\A\ 2 - \H 2 ) + i(l - u;)# 2 + ifojv 
< -arCl+^^e^ + OCe- 3 */ 2 ), 

where c > 0. Since the right-hand side of (|4.20p falls-off like t~ a , it follows from a standard 
maximum principle argument that < C(l + i) -cr where C depends only on the ambient 

geometry, and hence < Ce*/ 2 . Moreover since all the coefficients of (|4.20p and their 

derivatives are uniformly bounded, we may apply the local Schauder estimates on domains f2*r = 
dM x (T - T, T + T) to find 

II ^ llc 3 + a > 3 / 2 + a / 2 (n T/2 ) 

Ct (|| (1 + ty^e-^H- 1 ^ + x ~ w- 1 ) || cl+a , 1/2+a/2(nr) + || z || c °(n T )) , 
|| / \\ C k+ a ,k/2+ a /2 {nT) = sup F^/l + max [V r 5 t s /] a , a / 2 

n T r +2s<k r+2S ~ k 

[f]a,/3= SUp [f{-,t)] a ,dM + SUp [f(x,-)]p( To -T,T +T) 

T -T<t<T +T x&dM 

are Holder norms with exponent a in the space variable and f3 in the time variable. The desired 
bounds (14.16P now follow from the Schauder estimates, in light of the C° estimate that has already 
been established. 

■ 

As mentioned above our choice for z will ensure that, after perturbation, strict inequality holds 
in (|4.4p and (|4.5p . Choose X to solve 

(4.21) Ax = -2e 2 e~ 2t on M, x\dM = 0, x -> as |xrr| — > do, 

and choose to sufficiently small (depending only on the ambient geometry) so that equality cannot 
occur in (|4.4p or (|4.5p for < t < t . It follows that if Aw^ 1 = -e 2 e~ 2t then 

Aw' 1 = Aw' 1 + n(A x - Aw' 1 ) = -e 2 e~ 2t + fi(e 2 - 2e 2 )e~ 2t < -e 2 e~ 2t . 

Of course outside the set on which Aw -1 = — e 2 e~ 2t we may simply choose \i sufficiently small in 
order to obtain the same conclusion. Next observe that it is not possible for w~ x = 1 since w -1 
is super harmonic; (in fact this is the reason for including the superharmonic condition). Thus we 
need only consider the case when w -1 = 1 + ee~*. If this happens then 

vj- 1 = w- 1 + n(x - w- 1 ) = 1 + ee~ l - fi(l + 0(e)) < 1 + ee~\ 

assuming that e is small. Lastly consider the derivative estimates in (I4.4p . Pick local coordinates 
and let f3 be a multi-index of order I, so that the components of V w -1 are V^w- 1 . Let dp = ±1 
be such that IV 3 ui -1 ! = JgV^w -1 , and similarly let 5p = ±1 be such that [V^u^ 1 ! = dpS/^w" 1 . 
Notice that if 

sgnV^w" 1 ^ sgnV^uT 1 = sgn[(l - n)V p w~ x + 0(6,^V 2 *)] 

then 

[V^w^l < 2b l fie 2 e- 2t 
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< 

where 
and 



for some constant b\. Now suppose that {V^w 1 \ = 5aiee 1 with 5 > a, 1 ^ 2 . Then we have 



i 



= Sense' 1 + (fy - tf/OV^uT 1 - ^{Saiee' 1 + (8p - Sp)V p vT x + 0(bie 2 e~ 2t )) 
= 5aiee~ l - /j,(5a l £e~ t + 0(fye 2 e -2 *)) 
< (1 - n/2)SaiEe~ % 

— 1/2 

if e is small and the sequence a/ is chosen appropriately. If 5 < a l then the same argument 
shows that | V^w — 1 1 < ^a^ee - * + 0(bi/j,e 2 e~ 2t ). Thus strict inequality in the derivative estimates 
of (|4.4p is achieved by taking e small and a/ large. This completes the perturbation argument for 

(ED- 

We now consider (|4.5p . In order to calculate the first variation observe that (|4.17p and (|4.18p 
imply 

Dw~ x [z] := lim ^(uT 1 - w~ x ) = Hd t z - w^HDH' 1 ^] = X ~ uT 1 , 
/U— »0 

so that 

Dlog(u;fO[a!] = -{wH){H- x DuT x [z\ + w^DH^z]) = -wHd t z. 

Note also that 

Dt[z] = zHd t t = zH. 

We then have 

D(\\7\og(wH)\ 2 + eK - w' 1 + 1 - eee^)[z} 
= -2A(\7\og(wH),\7log(wH))z - 2(Vlog(wH),V{wHd t z)) 
+eDK[z\ + w' 1 - x + OezHe- 1 . 
In the case of equality in (|4.5p . the estimates for w , H, \A\, and z show that 
\A(\7log(wH),\7log{wH))z\ < |^||Vlog(wif)| 2 |z| 



< C(w~ l - 1 + e\K\ + Oee 



— i -r - I'M T " . <; *) 

i 



and 

|(Vlog(wff),V(wiIdt,z))| < |Vlog(wtf)|u;tf(|Vd t z| + |d t .zj|V log(wtf )|) 

< C^ -1 - 1 + e\K\ + Oee- 1 ) 1 ' 2 

< Ce 1 '^' 2 . 

Moreover the variation of the Gauss curvature is given by 

DK[z] = Vij(zA ij ) - A(zH) - zHK, 
and therefore in a similar manner 

el-D-KM + 9£\z\He~ t < Ce. 

Lastly since w~ x — x > 1/2 it follows that for small e the first variation is positive. This completes 
the perturbation argument for (|4.5p and thus completes the proof. 



Remark 4.3. Theorem 4.1 implies that any static vacuum solution (M, g, u) G £s whose boundary 
dM satisfies the non-degeneracy condition (jl.6p is outer-minimizing in M, at least if (M,g,u) is 
in the component of the standard solution. Note that the non-degeneracy condition is local at 
dM; it depends only on the 1 st order behavior of the geometry at dM, while the outer-minimizing 
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property is global in M. To our knowledge, this criterion for outer-minimizing surfaces is new even 
for the case of the simplest static vacuum solution, i.e. M 3 . 

Of course, there are many boundaries which are not outer-minimizing. As a simple illustrative 
example, let T 2 be a torus of revolution embedded in M 3 with H > 0. One may remove a (small) 
essential annulus from T 2 and smoothly attach two embedded discs to obtain a 2-sphere S 2 with 
H > 0. This construction may be performed to obtain a curve {S 2 )t, t € [0,1], of positive mean 
curvature spheres which for t < ^ are embedded and for t > ^ are immersed, with a single self- 
intersection point of the discs at t = |. For t < ^ but close to ^, the embedded spheres (S 2 )t are 
not outer- minimizing in R 3 , (and don't satisfy the non-degeneracy condition (11. 6|) ). 

This passage from embedded to immersed behavior also shows that the full boundary map Hb 
on £ + is not proper. It is easy to see that the embedded spheres {S 2 )t give solutions (M,gfi at , 1) 
which are in the component £ + containing the standard round solution. 

We mention again that it is unknown if there are other components (£ nd )' of Il^ 1 (B nd ), (or 
other components {B nd )' of boundary data in Met m,a (dM) x C™~ 1,a (dM) satisfying (jl.6p ). not 
containing the standard solution. Theorem 4.1 implies that either all solutions {M,g,u) in such 
components have a global positive mean curvature foliation close to IMCF, or none do. 

Remark 4.4. The proof of Theorem 4.1 shows that the lower bound (14.161) on the mean curvature, 
as well as the corresponding upper bound, is independent of e, and depends only on the ambient 
geometry of (M, g, u). Hence one may let e — > and it follows that (M, g, u) £ £ nd admits a global 
solution to the (pure) IMCF (|4.3p . with w = 1. We have not been able to prove this directly, 
since one does not have a direct analog of the non-degeneracy condition (|4.5p for the IMCF itself. 
Instead, one may view the flow (|4.4p as a regularization of the IMCF. 

5. Curvature estimates and Properness of IL3. 

In this section, we continue to work on (M,g,u) in place of (N,g^). Let injoM denote the 
injectivity radius of the normal exponential map from dM in M. The main result of this section is 
the following. 

Theorem 5.1. For (M,g,u) E £ nd , one has a global pointwise estimate 

(5.1) |-R|<A, 

on M, where A depends only on bounds for the Bartnik data (7,-ff) € B nd , m > 2. Moreover, at 
dM, one has the bounds 

(5.2) \A\<A, inj dM >A-\ 



The estimates (15. lj) . (]5.2p also hold for higher derivatives of R and A, up to order m — 2, m — 1 
respectively. 

Proof: For points in the interior of M, of bounded distance away from dM, this follows directly 
from the apriori interior estimates in [Anl| which state 

(5.3) \R\(x)<*— t \dlogu\(x) < A 



t 2 { x y 1 6 lv ; - t{xy 

where t{x) = dist(x, dM), where A is an absolute constant. Similar (scale-invariant) estimates hold 
for all higher derivatives of R and logn. So one only needs to consider the behavior near dM. At 
dM, the Gauss and Gauss-Codazzi (constraint) equations are given by: 

(5.4) \A\ 2 - H 2 + s 7 = -2R NN , 



(5.5) 



5(A - Hj) = -u~ 1 D 2 u(N, •). 
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Also, -.Raw = —Ric(N, N) = -u _1 AA(u) = u _1 (A 9 m« + HN(u)), so that 
(5.6) u{\A\ 2 - H 2 + s 7 ) = 2(A 9M u + HN(u)). 

From (|5.4p . a bound on at dM gives immediately a bound on \A\ at <9M, given control of 
(7, if). Similarly, a bound on \R\ on M gives a lower bound on the distance d con to the conjugacy 
locus of the normal exponential map expg m ■ 

Now, again under a bound on \R\, the outer-minimizing property (|4.ip from Theorem 4.1 implies 
a lower bound on the distance 5qm to the cut locus of expgM- To see this, suppose that 5qm « 1 
but A in (15. ip is bounded, A ~ 1. Then since d con is bounded below, there is a geodesic 7 of 
length 5qm in -W meeting dM orthogonally at points pi,P2- Let T be the boundary of the tubular 
neighborhood of 7 of radius r. Then T intersects dM in the boundary of two discs D\, D2 of radius 
approximately r, (for r small). If 5qm « r, then areaT < area{D\ U D2). Further, T C M is 
homologous to Lq U Z?2 in M. Removing then D\ U -D2 from dM and attaching T shows that dM 
is not outer-minimizing in M, giving a contradiction; (compare with Remark 4.3). 

Thus it suffices to obtain a curvature bound at or arbitrarily near dM. The higher derivative 
estimates may then be obtained by standard elliptic regularity methods. The proof of (15. ip is by a 
blow-up argument. If the curvature bound in (]5.ip is false, then there is a sequence (M, gi,Ui,Xi) G 
£ nd with bounded Bartnik boundary data such that 

\R 9i \(xi) -> 00. 

Without loss of generality, assume that the curvature of gi is maximal at X{. We then rescale the 
metrics gi to g[ so that \R\ is bounded, and equals 1 at Xi, 

(5-7) \R g ,\( Xl ) = l, \Rg,\( yi )<l, 

for any yi € (M,g'-). One may also need to rescale the potential it. For reasons that will be clearer 
below, choose points yi G M such that dist g '(yi,dM) = 1 and dist g <_(yi, xi) = 1, and rescale Uj so 

(5.8) Uifc/i) = 1. 

The sequence {M^g^Ui) has uniformly bounded curvature and uniform control of the boundary 
geometry, (boundary metric, 2 nd fundamental form and normal exponential map). By (|5.8|) and the 
Harnack inequality, the potential U{ is also uniformly bounded in compact sets. It follows from the 
convergence theorem in [ATJ for manifolds-with-boundary that a subsequence converges weakly, 
(i.e. in C 1,a ), to a C ' a static limit (X,g,u,x) with boundary (dX, 7,u). By the normalization in 
(|5.7p . the limit (A, g) is complete, (without singularities). Since dM is outer-minimizing in (M,gi), 
the C° convergence to the limit implies that d X is outer-minimizing in A: if D is any compact 
smooth domain in d X and D' is a surface in A with dD' = dD, then 

(5.9) areaD' > areaD. 

One has <9A = M 2 , the boundary metric 7 is flat, ff = 0, so dX is a minimal surface in A. One 
has u > in the interior of A, (by the maximum principle), but may have u = somewhere or 
everywhere on dX. The bound (|5.7p and the static equations imply that u« is uniformly bounded 
up to dX, within bounded distance to xi and the limit potential u extends at least C 1,a up to dX. 

We will prove below that the convergence to the limit is smooth, so that in particular 

(5.10) \R\(x) = l, 

where x = lim x j and R = Rx ■ 

On the blow-up limit (A, g), (|5.6p holds and becomes 

\u\A\ 2 = A dM u, 
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on dX. This equation holds weakly on dX with u G C l,a {dX); elliptic regularity then implies it 
holds strongly, and u € C 3,a (dX). Since u is harmonic, u is thus C 3 ' a up to dX. Also, by the 
Riccati equation N(H) = -\A\ 2 - R NN = -\A\ 2 + ^{\A\ 2 - H 2 + s 7 ), so that 

(5.11) N{H ) = ^\A\ 2 + H 2 - Sl ). 

This holds pointwise on the blow-up sequence (M, g[,Ui) and since s 7 — > and -ff — > for it 
follows that N(H) is defined pointwise on the limit dX and on <9X, 

(5.12) N(H) < 0, 

with equality on any domain only when ^4 = 0. 

Since dX is minimal, (|5.12p and the outer-minimizing property (|5.9p imply that 

(5.13) N(H) = 0, 

on dX. In more detail, (15. 9p and the fact that i7 = on 9 A implies the 2 nd order stability of dX, 
in that the 2 nd variation of the area of dX is non-negative. Thus, for all / of compact support on 
dX, one has 

(5.14) / (|d/| 2 + / 2 iV(i/))>0. 

Jdx 

Choose / = f R>s (r) such that / = 1 on D(R) C dX = R 2 and, for r > R, f = (log r-log 5) /(log R- 
logS), for S » R » 1. One may choose i? and 5 sufficiently large such that f dx \df\ 2 < e, for 
any given e > 0. This together with (|5,12p implies (|5,13p . 

It follows that A = and hence by the Liouville theorem on I 2 , u = const on dX. Using 
the divergence constraint (|5.5p . we also now have = 5(A — H^) = —u~ 1 D 2 (N, •), and so = 
D 2 (N, •) = dN(u) - A(du) = dN(u), so that N(u) = const. 

Thus, the full Cauchy data (7, u, A, N(u)) for the static vacuum equations is fixed and trivial: 7 
is the flat metric, A = and u, N(u) are constant. Observe that this data is realized by the family 
of flat metrics on (M 3 ) + with either u = const or u equal to an affine function on (IR 3 ) + . 

Suppose u = const > on dX. The static vacuum equations are then non-degenerate up 
to dX. The unique continuation property for Einstein metrics with boundary, cf. [AH], implies 
that the Cauchy data uniquely determine the solution locally. Alternately, since the static vacuum 
equations are non-degenerate up to dX and since the boundary data (dX, 7, H ) are real-analytic, 
elliptic regularity implies that the solution (M, g, u) is real-analytic up to dM. Such solutions are 
uniquely determined (locally) by their Cauchy data. Hence, the limit (X,g,u) is flat in this case. 

Moreover, the convergence to the limit is smooth everywhere. This again follows from non- 
degeneracy and ellipticity. Briefly, the potential equation Ait = gives a boost on the regularity of 
u, (given background regularity on g). One substitutes this into the main static vacuum equation 
uRic = D 2 u, giving thus a boost to the regularity of Ric, inducing then a boost to the regularity 
of g. This in turn further boosts the regularity of u via the potential equation. Bootstrapping 
gives C m ' a convergence, up to the boundary, in regions where u > 0, (given that u is C m ' a at the 
boundary). In sum, one has a contradiction to (|5.10p . 

Thus, suppose instead 

(5.15) « = on dX. 

This situation is more complicated. It is also more difficult to prove smooth convergence in this 
situation (one may have x in (|5.10p at dX). Moreover, there are in fact non-flat static vacuum 
solutions with flat Cauchy data as above with u = 0on dX, (so-called toroidal black holes, cf. [P] 
and pp). Thus the unique continuation results used above are false in this degenerate situation 
where the boundary becomes characteristic. 
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We observe first that the solution (X, g) is still real-analytic up to dX in this case. This follows 
since the 4-metric g 4 = u 2 d9 2 + gx is Einstein (Ric = 0) and is C 1,a up to the horizon or vanishing 
locus dX = {u = 0}. Elliptic regularity for the Einstein equations then implies that g is real- 
analytic, and hence so are u, gx up to dX. 

To see this in more detail, let U be a chart neighborhood of dX in X, so that U is diffeomorphic 
to a half-ball in IR 3 with boundary a disc D 2 C M 2 . Over each p £ U, one has a circle of length 
2ttu(p), with u — > as p — > <9X. From the work above, it is C 1,Q up to <9X with iV(it) = const 
at <9X. Note that TV (it) 7^ at <9X. For if TV (it) = at <9X, since also u = at OX and it is 
harmonic (Au = 0), the unique continuation property for harmonic functions implies that u = in 
X, giving a contradiction. By rescaling u if necessary, one may thus assume that N(u) = 1 at dX. 
This implies that the 4-metric g 4 , defined on S 4 \ D 2 extends to a C 1,a metric on the 4-ball B 4 . 
(The coordinate 9 is an angular variable in M. 2 in polar coordinates, shrinking down to the origin on 
approach to dX). It is well-known, cf. [Bej . that any C 1,a weak solution to the Einstein equations 
is real-analytic (in harmonic or geodesic normal coordinates), which gives the claim above. 

To prove the limit is in fact flat, and that one has strong convergence, we need to use the outer- 
minimizing property again. Thus, first note that (15. lip holds everywhere on the limit (X,g) near 
dX, not just at dX; here A is the 2 nd fundamental form of the level sets S(t) of t = dist(dX, •), 
etc. We have already established N(H) = at dX, via the outer-minimizing property and the 
corresponding stability of the 2 nd variation operator (|5.14j) . Taking then the derivative of (|5.1ip in 
the normal direction gives, 

(5.16) NN(H) = -(A', A) - (A 2 , A) - H'H + \s' v 

where A' = Vat A. At dX, the first three terms vanish while (s')k = —A(trk) + 55k — (Ric^, k) = 0, 
since k = 2 A = 0. Thus 

(5.17) NN{H) = 0, 

at dX and it follows that the 3 rd variation of the area of dX in the unit normal direction vanishes. 

Now choose / = f^s as following (15.141) with R, S large. Let Sff = expgx(tf(x)), where expsx 
is the normal exponential map of dX into X. Letting v(t) = areaStf, one has 

(5.18) v{t) = v(0) + \v"($)t 2 + ±1/" (0)i 3 + ±v""(0)t A + 0(t 5 ). 

The expansion (|5.18p is valid for all t sufficiently small, \t\ < 5q, with 5q independent of R, S, 
since the area and its derivatives are integrals of local expressions, and the local geometry of X is 
uniformly bounded in a tubular neighborhood of radius 1 about dX. By the 2 nd variational formula 
(|5.14p and (I5.13p . for any given e > 0, one has 

«"(0) < e, 

for R, S sufficiently large. For the same reasons via (|5.17p . 

v"'(0) < e. 

It follows then from the outer-minimizing property (|5.9p and (|5.18p that for R, S sufficiently large, 
one must have 

(5.19) v (0) > -s, 

again for any e = e(R, S) > 0. Using the vanishing of the lower order terms, one computes that 
(EOSD gives 

(5.20) / fNNN(H) - 6f 2 (df ■ df, A') > -e. 

Jdx 
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On the other hand, taking the normal derivative of (|5.16|) gives 

(5.21) NNN(H) = -(A', A'} - (H 1 ) 2 + ±4 

We have H' = N(H) = and (s' 7 ) 2A = 2AH + 255 A - (Ric^,2A). For s'^, one has (AH)' = 
A'H + AH' = and (Ric 7 , A)' = ((RiGy)', A) + {(Rioy), A') = 0. So at dX, = 55 A' . It follows 
then from (pT20l) - ([5^T| ) that for / = f R>s as above 

(5.22) / -f\A'\ 2 + f55A'-6f(df-df,A')>-e. 

JdX 

Integrating the second term by parts gives / (D 2 f 4 , A') = J (Af 3 D 2 f + 12 f 2 (df ■ df, A'). Using the 
Cauchy-Schwarz and Young inequalities, (|5.22p then implies, for any \x small, 

/ /VI 2 < /J /VT + cW f 2 \D 2 f\ 2 + C^[ \df\* + e. 

JdX JdX JdX JdX 

Choosing fj, small, the first term on the right may be absorbed into the left, while simple computation 
shows that the last two terms become arbitrarily small for R and S sufficiently large. It follows 
that 

(5.23) A' = 
and so NNN(H) = at dX. The Riccati equation 

(5.24) A' = V N A = -A 2 - R N , 

where Rn(V, W) = (R g (V, N)N,W), thus gives i?jv = at dX, and so via the Gauss and Gauss- 
Codazzi equations R g = at dX. Thus the full ambient curvature vanishes at dX. 

One can now continue inductively in the same way to see that A and R vanish to infinite order at 
dX. A simpler method proceeds as follows. The Riccati equation (15. 24ft holds along the level sets 
S(t) of t = dist(dX, ■). Since s g = and dimX = 3, Rjy = — * Ric, i.e. Rn(v,v) = —Ric(w,w), 
where (N,v,w) are an orthonormal basis. Via the static vacuum equations, this gives Vat A = 
—A 2 + n _1 * (D 2 u). Rescale u if necessary so that N(u) = 1 at dX and set v = u — t. Since 
A = D 2 t, one then obtains 

(5.25) V N A = -A 2 + u _1 (*A) + u~ 1 (*D 2 v). 

This is a system of ODE's for A, singular at dX = M 2 , but with indicial root 1. From the work 
above, we have v = 0(t 2 ) and A = 0(t 2 ). Writing A = t 2 B and substituting in (15.250 shows that 
A = 0(t 3 ). Also, by the computation following (I5.21|) . s 7 = 0(t 3 ) on S(t), and hence using the 
scalar constraint (|5.4p and the relation Rnn = u~ 1 NN(u), this in turn implies v = 0(t 3 ), and 
so on. It follows that (g, u) agree with a flat solution to infinite order at dX. Since the solution 
(X,g,u) is analytic up to dX, it follows that (X,g,u) is flat, as claimed. 

Next we claim that one has strong convergence to the limit, so that (|5.7j) is preserved in the 
limit, i.e. (15. 10H holds, contradicting the fact that the limit is flat. Note first that if x in (I5.10P is 
in the interior of X, then strong (C°°) convergence is immediate, by the interior estimates (|5.3p . 
i.e. their higher derivative analogs. Thus we may assume that x £ dX. 

From the work above, we know that \R\ is uniformly bounded everywhere on (M, g^) and \R\ — > 
everywhere away from dM — > dX, so \R\ jumps quickly from 1 to near X{. The main point is to 
prove that 

(5.26) \R\(x)^0, 

for all x £ dM — > dX; it is then easy to prove that \R\ — > on (M, g^), cf. (|5.35p below. To prove 
(I5.26p . note that the estimates (I5.1ip - (I5.23|) above at dX also hold on the blow-up sequence at dM, 
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(since H — > and s 7 -> on dM). It follows then by these arguments that for any R < oo and 
D(R) the i?-ball about any base point y, t G (dM, 7^) converging to y G dJT, 

(5.27) / |^| 2 + |Viv^| 2 -> 0. 

To proceed further, consider the divergence constraint 5A = —dH — Ric(N, •) on (dM, 7-), as in 
(|5.5|) . The equations 5 A = xi) di^-A = X2 form an elliptic first order system in 2-dimensions, and 
so one has elliptic estimates. Since H -> G C™" 1 -" on (<9M,7^), dfT -> in C" 71 " 2 ' . Also, 
by assumption (|5.7|) . Ric(N,-) is bounded in L°°. It follows then from elliptic regularity that A 
is bounded in L 1,p . By the scalar constraint (|5.4p . -Ratat is then also bounded in L 1,p and since 
trR N = R NN , it follows from (I5T271) and (|5T24j) that 

(5.28) ,4^0 and R NN -» in Cf oc (dM). 

Next the Einstein equation Ric gN = on (N,gjy) implies that 5nR 9n = 0. Hence 

= (5 ffJv i?)((-, iV)AT, •) = 5 gjv (i?(-, iV)AT, •) + 2i?(e a , V ea iV)iV = 5 aN (R(-, N)N, •), 

since one may choose a basis in which V ea iV = A(e a ) = X a e a . Let V be the unit vertical vec- 
tor, and note that V V V = -dv, where v = logu. Then V V (R(V, N)N, ■) = R(dv,N)N, while 
VAr(i?(A r , N)N, ■) = 0. Hence, for Rn as in (|5.24p . these computations on dM give 

5Rn = —RN(dv), 

where the divergence 5 and Rn are taken on (dM, 7-). Since i? is bounded in L°° and fri?Ar = Rnn 
is bounded in L 1 ^, elliptic regularity gives 

(5.29) \\R N \\ L x, P <C\\dv\\ LP , 

again on compact domains in dM converging to a compact domain in dX. 

To control dv in (|5.29p . recall that by (15. 7p the ambient curvature R is bounded, and so R 
restricted to dM is also bounded. Via the static vacuum equations (jl.ip . this implies that 

(5.30) N(v)A + u~ l D 2 u 

is bounded on (dM, 7-), where D 2 u is the Hessian of u\qm '■ dM — » M + . Now we claim that each 
term in (|5.30p is bounded, i.e. there exists K such that 

(5.31) {u^D 2 ^ < K, \N(u)A\ < K, 

pointwise, on domains converging to a bounded domain in dX. To prove (|5.3ip . suppose instead 
that |u _1 L» 2 u| — > 00 at some sequence of base points yi — > y G Without loss of generality, we 
may assume the points yj realize the maximum of \u~ 1 D 2 u\ on D yi (W) C (9M, 7^), (possibly up to 
a factor of 2). As before, one may then rescale the metrics g[ further to g'( so that \u~ l D 2 u\(yi) = 1 
and hence the full curvature R — > in this scale. Also, renormalize u if necessary so that u(yi) = 1. 
Note that \du\(yi) must be bounded. For if \du\(yi) were too large, it follows, (e.g. by a still further 
rescaling), that u would be close to an affine function on M 2 and hence u would assume negative 
values in bounded distance to Since u > everywhere, this is impossible. Thus, by integration 
along paths, u is bounded in L 1,oa in the scale g'(, within bounded distance to y-i. 

Now working in the scale and normalization above, from divergence constraint f)5.5|) . one has 
— u5(A — Hj) = dN(u) — A(du). Since A is bounded in L > p and u and du are bounded in L°° , it 
follows that dN(u) is bounded in LP . Thus, N(u) = c + 4>, where (j) is bounded in L ,p . Here c is a 
constant which may, and in fact does, go to 00, in the n-normalization u(yi) = 1 above. 

The trace equation (|5,6p in this scale and normalization gives 



(5.32) 



Au + H(c + 4>) =uf, 
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where / is bounded in L 1,p . Also, He is bounded, since N(u)A is bounded via (|5.30p . and so 
He — > d, for some constant d on dX. Since <j> is also bounded in L 1,p , it follows that u is bounded 
in L 3 ' p , and hence (in a subsequence), u converges in C 2 ' a to its limit on dX = M?. Hence D 2 u 
converges in C a to its limit on M?. On the limit, since H — > 0, the trace equation (|5.32j) becomes 

(5.33) Au + c' = 0. 

If d = then since u > 0, u = const and hence .D 2 it = 0, giving a contradiction. If d ^ 0, then 
again since u > 0, one must have c' < and it is a quadratic polynomial on R 2 . Since u is harmonic 
on (M 3 ) + and d < implies JV(it) = const < 0, this is also easily seen to be inconsistent with the 
requirement u > everywhere. This proves (|5.3ip holds. 

Returning to (|5.29|) . we may work in the normalization above that u(yi) = 1 and then (|5.3ip 
implies that dv = u~ x du is bounded in L°° in bounded domains about Hence by (|5.29j) . Rn is 
bounded in L 1,p and so bounded in C a . Since Rn — ► in L 2 locally on dM , one has 

(5.34) R N -> in Cg c (dM). 

This is the main part of the estimate (|5.26p . 

Next, one needs the same result for the Ric(N,T) term. To do this, take the normal derivative 
of the scalar constraint (|5.6p . to obtain 

AN(u) + A'u + N(H)N(u) + HNN(u) = ±N(u)[\A\ 2 - H 2 + s 7 ] + \uN[\A\ 2 - H 2 + s 7 ]. 

Since 7' = 2^4, a standard formula for the variation of the Laplacian, (cf. |Be| ). gives \A'u = 
-(D 2 u,A) + (du,(3(A)) which is bounded in L°°. Also the terms N(H), H and NN(u) = -f [|^| 2 - 
i? 2 + s 7 ] all go to in L°°. Moreover, via (f5T34"|) above, uN[\A\ 2 - H 2 + s^} is bounded in L°°. Also, 
as in the proof of (|5.3ip . N(u) = c + with bounded in L 1,p . Thus it follows from the above by 
elliptic regularity that dN(u) is bounded in L 2,p and so converges in C 1 '" on dM. Since 

-uRic(N, ■) = dN(u) - A(du), 

it now follows that Ric(N, ■) converges to its limit, necessarily 0, in C a on dM. Lastly, 5 A = 
—dH — Ric(N, ■) — ► in C a and hence by elliptic regularity, A — > in C 1 '" so that via the Gauss- 
Codazzi equations dA(X, Y, Z) = (R(N, X)Y, Z) — > in C a ; here ci is the exterior derivative. 
Combining the computations above now proves the estimate (|5.26p . 

To complete the proof, we have \R\ ^ on dM -> dX. On the 4-manifold N = M x u S 1 , the 
Einstein equations give the inequality 

(5.35) AatI^I + cI^I 2 >0, 

where c is a fixed numerical constant, A^r is the 4-Laplacian and R 9N is the curvature tensor on 
N. One has \R gN \ = \R\, (up to a constant). We have proved above that R — > in L 2 locally 
on (N, (<7iv)i) and R — > pointwise at 9A^. It follows from the deGiorgi-Nash-Moser estimates 
for domains with boundary, (cf. [GT, Thm. 8.25]), that R — > pointwise on and hence on M, 
contradicting (|5.7p . This completes the proof. 

■ 

Next, we show that the potential function u is also controlled by the boundary data of n#. 

Corollary 5.2. For (M,g,u) £ S nd , there is a constant Uq depending only on the Bartnik data 
(7, H) G B nd such that 

(5.36) u < U , 

on M. Moreover, if H > Hq > on dM, then there exists Uq as above depending in addition only 
on Hq, such that on M, 

(5.37) u > Uq 1 . 
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Proof: Let 5(s) = {x G (M, g) : dist(x, dM) = s} be the geodesic 'sphere' about dM. Choose 
a fixed base point xo G 5(1) and suppose one has the bound 

(5.38) Cq 1 < u(xq) < cq. 

By Theorem 5.1, the geometry of the annular region A(^,2) about 5(1) is uniformly controlled by 
the Bartnik data (7, H) and so by integration of the static vacuum equations uRic = D 2 u along 
paths in A(^, 2), one has 

(5.39) Cq 1 <u<C , 

in A(j, |), where Cq depends only on cq and (7, if). 

To remove the dependence of Cq in (|5.39p on cq in (|5.38p . we need better control on the large-scale 
behavior of u. To do this, it is proved in [An2, Lemma 3.6], that there is a constant K, depending 
only on Co, such that for all s > 1, 

(5.40) sup \du\ < K^Vcis))- 1 , 

S c (s) 

where v c (s) = areaS c (s) and S c (s) is any component of the geodesic sphere S(s). We point out 
that (|5,40p holds for general static vacuum solutions, not only those in 8 nd for instance. The 
estimate ()5.40p is proved by studying the behavior of the harmonic potential log u on the Ricci-flat 
4-manifold (N,gw, log u) and then reducing to (M, gM,u). 
Consider now the conformally equivalent metric 

(5.41) g = u 2 g. 

It is well-known that the static vacuum Einstein equations (jl.ip are equivalent to the equations 
Ric = 2(du) 2 > 0, AgV = 0, v = logn. The metric g thus has non-negative Ricci curvature with 
harmonic potential v. These are exactly the properties used to prove (|5.40p . and a brief examination 
of its proof shows that (|5.40p also holds with respect to g, i.e. 

(5.42) sup \dv\g < K(v c (s))-\ 

S c (s) 

again with K = K(Cq). Since (M, g) is asymptotically flat and u — ► const at infinity, the area 
growth of geodesic spheres v(s) in (M, g) satisfies v{s)/s 2 — * u>2, where u>2 = areaS 2 (l). It follows 
then from the volume comparison theorem for Ricci curvature, (cf. |Pej ) . that 

(5.43) v(s) > uj 2 s 2 , 

for all s > 1. As above, by integration of (|5.42p along a geodesic ray starting from a suitable base 
point x\ G 5(1) out to infinity, one sees that (|5.36p holds then globally on M\B(1), with Uq again 
depending only on cq in (|5.38|) . Using the static vacuum equations, the same integration along 
paths gives such a bound within B(l). Thus, we see that (|5.36p follows from (|5.38l) . 
To prove (|5.38p . suppose one has a static vacuum solution (M,g,u) with 

(5.44) u(x ) = e. 

Renormalize u to u = u/u(xq), so that u(xq) = 1 and u — > e~ 1 at infinity. Then (|5.38p holds, and 
hence so does (|5.42|) - (|5.43|) . Again by integration along geodesies starting at xq and diverging to 
infinity, it follows that 

u < Ui, 

where U\ depends only on the boundary data of IL3. This proves the lower bound in (|5.38j) . 
To prove the upper bound, suppose instead 

(5.45) u{xq) = e" 1 . 

29 



Then again we renormalize u to u as above, so that now u — > e at infinity. This does not directly 
give a lower bound on e via (|5.42j) - (|5.43j) as above. However, one may proceed as follows. First, it is 
well-known that static vacuum solutions come in "dual" pairs, in that if (M, g, u) is a static vacuum 
solution, then so is (M,g,u) with g = u 4 g, it = u~~ l , cf. [An2| for instance. Then (|5.42p - (|5.43p hold 
for (M,g,u) which as before by integration gives an upper bound u < Ui at infinity. Since near 
infinity, u ~ e _1 , this again gives a bound on This completes the proof of (|5.36p . 

To prove (|5.37p . note that (|5.38p has been proved above, and hence by the maximum principle 
and normalization u — ► 1 at infinity, (|5.37p holds in the exterior region M \B(T). Thus, one only 
needs to consider the behavior near dM. For this, suppose (M,g,u) is a static vacuum solution, 
C 2 ' a up to dM with u>0onM = MU dM. If H > H > on dM, we claim that necessarily 

u>0 on dM. 

For if u = at some point z G dM, then by (|5.6p . Aqmu + HN(u) = at z. Since = u(z) is 
a global minimum for u, one has Aq^u > and by the Hopf maxmimum principle, iV(u) > at 
z. This gives a contradiction if H > 0. The same arguments prove the existence of a lower bound 
(|5.37p by a contradiction argument, taking a sequence and passing to a limit, using Theorem 5.1. 
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The previous results now lead quite easily to the following main result of this section. 
Corollary 5.3. The boundary map 

U B : 8 nd -» B nd , 

is proper. 

Proof: Let (M, gi,U{) be a sequence of static vacuum solutions in £ , with ile^i, «j) = (7$, 
Supposing (74, flj) — > (7, if) in B nd , we need to prove that the sequence (gi, Ui) has a subsequence 
converging in C m,a (M), modulo diffeomorphisms, to a limit (M,g,u) G £' nd . 

The curvature bound (|5.ip and control of the intrinsic and extrinsic geometries of the boundary 
metrics first implies the metrics gi cannot collapse within bounded distance to dM, i.e. there is a 
fixed constant iq > such that the injectivity radius of (M,gi) satisfies 

inj 9i {x) > i , 

for dist gi (x,dM) < K. By the compactness theorem in, for instance [AT] , it follows that a sub- 
sequence of (M,gi) converges in C m ' a , (and C°° in the interior), uniformly on bounded domains 
containing dM, to a limit (M',g). One has dM' = dM and g is a complete Riemannian metric on 
M', C m ' a up to dM and C°° in the interior. 

By Corollary 5.2, the potential functions u% also converge in C m ' a , (in a subsequence) to a limit 
potential function u on M', and the pair (g,u) gives a solution of the static vacuum Einstein 
equations. Since u = limuj, it follows that 

(5.46) u < U Q , 

on M', for C/o as in (|5.36p . Clearly the boundary metric and mean curvature of (M' , g, u) are given 
by the limit values (7, if). To prove that (M',g,u) G £ nd , one then needs to prove that (M',g,u) 
is asymptotically flat and M' is diffeomorphic to M. Note that since the convergence above is only 
uniform on compact sets, apriori there need not be any relation between the asymptotic structure 
of (M',g,u) and (M,gi,Ui) for any given i. 

The equation (I5.40P holds on each (M,gi,ui) and by Corollary 5.2, the constant K is uniform, 
independent of i. Moreover, as in the proof of Corollary 5.2, there is a geodesic ray a = Oi starting 
at any fixed base point in S(l) and diverging to infinity, such that on the component S c (s) of S(s) 
containing a, one has 

sup \dui\ < Ks~ 2 , 

Sc(s) 
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where K is independent of i. Since Ui is harmonic, by elliptic regularity, (and scaling), a similar 
estimate holds for higher derivatives of Ui, and via the static vacuum equations, it follows that 

sup \R gi \ < Cs~ 3 , 
S c (s) 

with again C independent of i. This means that the metrics (M,gi,Ui) become asymptotically flat 
at infinity uniformly, at a rate independent of i. For R sufficiently large, R > Rq independent 
of i, the geodesic spheres S(R) and annuli A(R, 2R) are close to Euclidean spheres and annuli, 
(when scaled by R -1 ), and hence the geometry is close to that of Euclidean space; there can be 
no branching or joining of different components of S(R) for R > Rq. This implies that the limit 
(M',g,u) has a single asymptotically flat end, and M' is diffeomorphic to M. 

m 

Remark 5.4. The results above also show that the boundary map Us is proper not only on the 
component £ nd but also its closure £ . In other words, if (M, gi,Ui) is a sequence of static vacuum 
solutions in £ nd with boundary data (7^, Hi), (Hi > 0, dHi ^ on {Ki < 0}), and (ji,Hi) — * (7, H) 
in Met m ' a (dM) x C m ^ 1 ' a (dM), then (M,gi,Ui) converges in C m ' a (in a subsequence) to a limit 
(M, g, u) in £ nd , with H > 0. 

To see this, note that the constant A in Theorem 5.1 does not depend on a positive lower bound 
on H, and so (|5.ip - (|5.2p hold for the sequence (M,gi,Ui) above. Of course we are using here 
the fact that dM is outer-minimizing on the sequence (M, gi,Ui). Similarly in Corollary 5.2, the 
upper bound Uq on Ui does not depend on a lower bound for H. One may then use the argument 
concerning (|5.44p to show that Ui cannot go to on M away from dM. The proof of Corollary 5.3 
also does not require a bound on H away from 0. 

It is also worth pointing out that the results of this section only require that dM is outer- 
minimizing in a neighborhood of arbitrarily small but fixed size about dM. 

6. Degree of rig. 

By |Smj . a smooth proper Fredholm map F : B\ — > B2 of Fredholm index between connected 
Banach manifolds B\, B2 has a well-defined degree (mod 2). Namely, if y is a regular value of F, 
then i ?_1 (y) is a finite set of points, and deg% 2 F is just the cardinality of F^ 1 ^!/) (mod 2). In fact 
if Bi and B2 are oriented, then F has a well-defined degree in Z. 

By Corollary 5.3, it thus follows that the boundary map 

U B : £ nd -» B nd , 

has a well-defined mod 2 degree, degz 2 TlB- We do not know if £ nd is orientable, and do not pursue 
this issue here. 

The main result of this section is the following: 

Theorem 6.1. For (M, dM) = (M 3 \ B 3 , S 2 ), one has 

(6.1) deg% 2 n B = 1- 

Proof: The proof is based on the black hole uniqueness theorem [I], [Rj, [BM| . that the 
Schwarzschild metrics 

(6-2) 9Sch(m) = (1 - — Y l dr 2 +r 2 g S 2 {1) , u = J 1 - — , 

r > 2m, are the unique AF static vacuum metrics with a smooth horizon 7i = {u = 0}. Of course 
the Schwarzschild metrics are not in £ nd , but instead lie at the boundary d£ nd . 

Consider for instance any sequence {(ft, Mi)} G £ nd for which iIs(ft,Uj) = (7^,^) — > (7,0) 
smoothly, with K 7 > 0. Clearly, {(g i: U{)} is a divergent sequence in £ nd . By Corollary 5.3 and 
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Remark 5.4, a subsequence of {(gi, Uj)} converges smoothly to a static vacuum limit (M, g, u). (Of 
course one may have « = 0on dM). On this limit, 

H = 0, 

at 9M, so that <9M is a minimal surface. From (|5.6p one has 2Aqmu = u(\A\ 2 + s 7 ) > 0, and hence 
it follows from the maximum principle that u = on dM. Via the static vacuum equations (|l.ip . 
this implies further that A = and iV(«) = const at <9M. The black hole uniqueness theorem then 
implies that any such limit is the Schwarzschild metric, and so unique up to scaling. Thus one has 
uniqueness for the boundary data (7,0), so that most all boundary metrics 7 cannot be realized 
with H = at dM, (the no-hair result). 
Given this background, suppose 

(6.3) degz 2 H B = 0. 

Then for any reg ular value (7, if) G B nd , the finite set U]^(-f,H), if non-empty, consists of at 
least two distinct static vacuum solutions (g , u 1 ), (g 2 ,u 2 ). The regular values of LTb are open and 
dense in the range space (by the Sard-Smale theorem). Choose then a sequence of regular values 
(7$, .Hi) — > (7+1,0) smoothly. Note that there exist such regular values and corresponding regular 
points, since there are boundary data (j,H) £ HB(£ nd ) arbitrarily close to the Schwarzschild data 
(7+1,0); for example, one can go in a bit from the Schwarzschild horizon. (We set m = 1/2 here). 

Let (gj,uj), (gf,u 2 ) be any pair of corresponding distinct sequences in LT^ 1 ^, Hi). By Corollary 
5.3 and Remark 5.4, the sequences (gj,uj), (gf,u 2 ) have C m,a convergent subsequences to limits 

(9loi u lo)i {gloria) m £ n and W the uniqueness above 

9L=glo= 9Sch{ni), 

with m = 1/2, with = u 2 ^ = u in (16. 2D . 

This implies that near gs c h, the boundary map IL3 is not locally 1-1, and so presumably DUb 
has a non-trivial kernel at gsch- (Note however that gs c h £ We claim this is impossible. To 
prove the claim, let 

9Sch = u 2 d9 2 + g Sc h, 
be the 4-dimensional Schwarzschild metric on M? x S 2 , and similarly let 

be the 4-dimensional static Ricci-flat metrics associated to (gf,vE). By Lemma 2.2, without loss 
of generality we may assume that each is in Bianchi gauge with respect to Qsch, so that, as in 
pT0]) - (l2TTT|) . 

for j = 1,2 and i sufficiently large. By the smoothness of the convergence above, one may write 
(6-4) g|= 0Sc , + e^ + O((4) 2 ), 

where L(k-) = and L is the linearized Einstein operator (|2.6p at Qs c h- The data gj, Qsch and k\ 
are all smooth, (up to the boundary). The forms k\ are only unique up to multiplicative constants, 
which will be determined by choosing so that the C 1,a norm of — Qs c h equals e\ . Thus the 
C 1,a norm of k\ is on the order of 1. Note that decays to at infinity, so it is basically supported 
within compact regions of M. Let e$ = max(e^,e?). Then 

er^fli ~Bi) = K t + 0(ei), 

where Ki = e~ 1 (e 2 k 2 — s\k\) — > k, where the convergence is in C 1,a ', (in a subsequence). As 
previously, we need to show that the convergence is strong, so that k / 0. This follows from a 
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standard linearization and bootstrap argument, as preceding (|5.15p . In more detail, dropping the 
index i, we have A g jU 3 = 0, so that 

A g 2 (u 1 — u 2 ) = (A g 2-A g i)u\ 

In local harmonic coordinates, the right side this equation is on the order of e in C > a , and hence 
by elliptic regularity, u\ — U2 is on the order of e in C 3,a . Substituting this in the difference of the 
static equations v?Ric g j = D 2 jU 3 and arguing in the same way shows that the difference g 1 — g 2 is 

then also on the order of e in C 3 ' Q . This proves the strong convergence. 
It follows that the limit form 

(6.5) K = (h,u'), 

is a non-zero C 1,a weak solution of the linearized static vacuum equations L(n) = at gsch an d 
since U B (g},uj) = U B (g 2 ,u 2 ), one has 

(6.6) ry> h = H' h = at dM, 

where 7^ = h T = Ii\qm- As discussed following (|5.15p . elliptic regularity implies that (h,u') is 
smooth and so in particular a strong solution. Below we will use the fact that the data (h, u') are 
in fact real-analytic up to dM, (again by elliptic regularity). 
We claim that 

(6.7) {h, v!) = on M, 

which will give a contradiction. This is of course a linearized version of the black hole uniqueness 
theorem. It is possible that (|6.7j) can be proved by linearizing one of the existing proofs of black 
hole uniqueness in [I], [R], [BMj . However, we have not succeeded in doing this and instead (|6.7p 
is proved in a manner similar to the proof of Theorem 5.1. 
First, the linearization of (|5.6p gives, at dM, 

u's 7 = 2Au'. 

Since s 7 > 0, the maximum principle implies that v! = at dM. Next we claim A' = 0. To 
see this, the vacuum equations give uRic = D 2 u and D 2 u = N(u)A + (D 2 u) T when evaluated on 
tangent vectors to dM. Taking then the variation and evaluating tangentially gives (uRic)' = 
so = (D 2 u)' = (D 2 )'u + D 2 u' . The first term on the right vanishes when evaluating tangentially 
and hence so does the second term. This implies = (N(u)A)' = N(u)A' + N{u')A = N(u)A' . 
Since N(u) = const ^ 0, it follows that A' = 0. Similarly, taking the variation of the divergence or 
vector constraint gives N(u') = const. 

Clearly N(u') = to', (up to constants). A simple examination of the proof of black hole unique- 
ness in [R] applied to an Einstein deformation as in (|6.5|) and satisfying (j6.6fl . shows easily that 
N(u') = at dM. (One does not obtain any further information, since the bulk data in the 
Robinson proof, via divergence identities, are quadratic in the deviation from Schwarzschild). 

Thus the variations (7' ,u' , A' , N(u')) of all the Cauchy data are trivial. As in the proof of 
Theorem 5.1, we use a bootstrap argument to prove that the data (h, u') vanish to infinite order at 
dM, in geodesic gauge. 

Thus, using geodesic normal coordinates near dM, write 

g = dt 2 +g t , 

where t(x) = dist(x,dM). We may assume, (by adding an infinitesimal deformation of the form 
5*V if necessary), that h preserves this gauge, so that ho a = 0, i.e. h{N, ■) = 0, N = dt, near dM. 
By the discussion above, we have u' = N(u') = at dM and similarly h = V^h = at dM, so 
that 

(6.8) u' = 0{t 2 ) and h = 0(t 2 ). 
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The variation of the potential equation Amu = gives 

(6.9) Au = -A'u = (D 2 u, h) - (/3(h), du), 

where (3 is the Bianchi operator, (cf. |BeJ). Since (3(h) = at dM, this gives Au' = at dM and 
hence NN(u') = at dM, so that 

(6.10) u' = 0(t 3 ). 

Next the linearization of the Riccati equation gives (V^A)' = — (A 2 )' — (Rn)' = — (i?jv)' at dM. 
One computes *Rn = —(Ric T ) = —u~ 1 (D 2 u) T , so (*Rn)' = u~ 2 u'D 2 u — u~ 1 (D 2 )'u — u~ 1 D 2 u', as 
a form on T(dM). It follows from (fBTTUj) that (*R N )' = at dN and hence (V N A)' = so that 

(6.11) h = 0(t 3 ). 
Next taking the normal derivative of (|6.9|) gives 

AN(u') = (V N D 2 u,h) + (D 2 u,V N h) + (X7 N (3(h),du) + ((3(h),V N du), 

which vanishes at dM and hence v! = 0(t 4 ). Substituting this in the linearized Riccati equation 
above and using previous estimates gives h = 0(t 4 ), and so on. It follows that (h,u ! ) vanish to 
infinite order at dM. Since (h,u') are real-analytic up to dM, (in geodesic gauge), this implies 
that h = u' = on M, i.e. (|6.7p holds. This completes the proof. 

■ 

Although the proof of Theorem 6.1 implies that DUg has trivial kernel at the Schwarzschild 
metric gsch, one does not expect this to be the case for the cokernel. In fact, one expects that 
CokerDUs is infinite dimensional, in that any boundary variation of the form (k, 0), where k is a 
variation of the boundary metric, is not tangent to a curve of static metrics with H = at dM. 
This amounts to the linearized version of the no-hair theorem, (which has not been proved as far 
as we are aware). In particular, we expect DUb is not Fredholm at gsch- 

Remark 6.2. We point out that ii£ nd is orientable, then the proof above shows that degjjls = ±1 ; 
depending on the choice of orientation of £ nd , (and the target space). 

Remark 6.3. It is interesting to compare Theorem 6.1 with the bounded domain case. Thus 
suppose (M,dM) = (B 3 ,S 2 ) and consider the space £ = £ m ' a of C m,a static vacuum solutions 
on M. This is again a smooth Banach manifold for which the boundary map 11^(9) = (~f,H) is 
Fredholm, of index 0. By general theory, Ub thus always has a local degree. 

We claim that the local degree of near the standard flat metric solution (g, u) = (gjiat, 1) on 
B 3 (l) is 0, 

degi oc U B = 0. 

Thus, consider static vacuum metrics (g, u) near the flat data, with induced boundary data (7, H). 
Since > 0, there is a unique isometric embedding (dM, 7) — > M 3 , with induced mean curvature 
Hq. By a basic result of Shi- Tarn [ST] . 

I HdV^ < [ H dV 1 , 
JdM JdM 

so that the functional J dM (H — Ho)dV^ on £ has a local maximum at (g, u) = (gfi a t,l). In 
particular, the boundary data (7+1, H), for any H with H > 2, is not in ImHB, so that n# is not 
locally surjective. Hence, the local degree is 0. 

It also follows that the solution (gfiat-, 1) ° n B 3 (l) is a critical point of Ub- In fact there are curves 
(g t ,ut) of static vacuum solutions with (50,^0) = (fl 'flat A) and U B (g t ) = (7+i> 2 ) = Kb ((g fiat A))- 
For this, one can just take g% = gfiat and Uf = 1 + tz, where z is any afiine function on M 3 . 
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Similarly, the main result in [Mi] also shows that the standard flat data (gfiat, 1) on M 3 \ S 3 (l), 
with boundary data (7+1,2) is a critical point of n^. Thus, the standard flat boundary data 
(7+1, 2) is not a regular value of Us, for either the interior or exterior boundary map lie. 

As noted in the Introduction, Theorem 6.1 implies Theorem 1.1, since any smooth proper Fred- 
holm map of non-zero degree is surjective. 

We conclude this work with several remarks and questions worthy of further study. First, it 
would be interesting to know when solutions (M, g) of the static vacuum Einstein equations (|1.1|) 
are uniquely determined, (up to isometry), by the Bartnik boundary data (7, H). This is unknown 
for solutions in 6 ; in fact it is even unknown for solutions close to the standard solution (IR 3 \ 
B 3 (l) 

,9fiatA), cf. Remark 6.3 above. 
Note that the mass m of the solutions in £ nd is not necessarly non-negative, but can take on all 
values in BL It is an interesting question to determine conditions on the data (7, H) for which the 
corresponding solution(s) {M,g) have m > or m > 0; cf. [H, Prop. 2.1] for some results on this 
issue. For static vacuum solutions, the mass m can be computed concretely at dM as 

m = -3- / N(u)dAj. 

4tT JdM 

Finally, recall that Bray [Brj has suggested replacing the "no horizon" condition in the definition 
of admissible extension in (|1.5jl by the property that dM is outer-minimizing in (M,g). Con- 
ceptually, this meshes well with the results here, (since all solutions in £ nd have outer-minimizing 
boundary). However, as discussed previously, it is not easy to determine from the Bartnik boundary 
data alone which solutions have outer-minimizing boundary. 

7. Appendix. 

Here we derive the first variation of energy formula (|4,7p - (|4.8p used to obtain the bounds for 
the mean curvature for the generalized IMCF in §4. The first task is to derive the appropriate 
evolution equation. Recall that the flow is given by dt = (wH)^ 1 N; we will use the notation of §4. 
Therefore 

d t H = -AiwH)- 1 -(wH)-\\A\ 2 + R NN ) 

= (wH)- 2 A(wH) - 2(wH)- 1 \Vlog(wH)\ 2 - (wH)- l (\A\ 2 + R NN ), 

and, since u is harmonic on (M,g), 

wH 2 dtu = —Au — Rnnu. 

Let P(u) be a positive function of the potential. Then 

(7.1) d t {wHP) = wP[(wH)- 2 A{wH) - 2(wi7)" 1 1 V log(u;if)| 2 

-(wH)-\\A\ 2 +R NN ) +w- l N(logw) - {wH)- 1 (log P)'(Au + R NN u)]. 

Now consider 

f-% [ ^(x)(f(t)wHP(u)G(u)) x 

= [ X^(wHP) x - 1 G x d t (wHP) + X(log f)'^(wHPG) x 

+ X^iwHPG^iwH^ilogGYuNilogu) + ^(wHPG^w- 1 , 
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where G(u) is another positive function of the potential, to be chosen below. According to (|7.ip 
there are two terms which should be integrated by parts, namely 

f if>(x)(wHP) x - 1 G x wP(wH)- 2 A(wH) = - [ i)w{PG) x {wH) x - 2 Q x , 

where 

Qi = (\-3)\Vlog(wH)\ 2 + (Vlog(wtlj),V\og(wH)) 

+A(logP)'u(Vlogu, Vlog(uPT)) + A (log G)'u(V log u, V log(ioiJ)), 

and 

-/ ^(x)(wHP) x - 1 G x PH-\\ogP)'Au= [ iJw(PG) x (wH) x - 2 Q 2 , 

where 

Q 2 = (A - l)(logP)''u(Vlog(u;i7), Vlogu) + ((A - l)P' 2 p- 2 + P'^^fVlogul 2 
— (log P)'u(V logu, V log(wH)) + (logP)'u(Vlogu, Vlogw) 
+A(log G)'(log P)'u 2 1 V log u\ 2 + (log P) r u(V log V, V log n) . 

We now have 

/ Hx)(f(t)wHP(u)G(u)) x = [ \4,w(PG) x (wH) x - 2 Q 3 , 

JT, t JT, t 

where 

Q 3 = -(A- l)|Vlog(u»P)| 2 - (2(logP)' + A(logG) / )'u(Vlogn,Vlog(^)) 

+((log P)" + A(log P)' 2 + A(log G)'(log P)')w 2 1 V log u\ 2 - \A\ 2 + X^H 2 

-(1 + (logP) / n)Pj ViV + w{\ogf)'H 2 + HN(\ogw) + fl"(log G)'uN(\ogu) 

-(Vlog(wip),(Vlog(wH) - (logP)'uVlogu)). 

In order that the quadratic form involving first derivatives of log(u>-fPj and logu be definite, we 
choose G = P~ 2 . Then upon completing the square the quadratic form becomes 

((logP)" - (logP)' 2 )n 2 |Vlogn[ 2 - (A - 1)|V log(wH) - (logP)'nVlogu| 2 . 

Lastly setting F = P~ l yields the desired expression (I4.7p - (|4.8p . 
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